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FOREWORD 


This  is  the  final  Technical  Report  prepared  by  the  Cal span 
Corporation.  The  effort  was  sponsored  by  the  Air  Force  Aero- 
Propulsion  Laboratory,  Air  Force  Systems  Cormand,  Wright-Pat terson 
AFB,  Ohio  under  Contract  F33615-73-C-2046  for  the  period  1 May  1973 
to  31  May  1976.  The  vsork  herein  was  acconplished  under  Project  3066, 
Task  306603,  Work  Unit  30660334,  "Investigation  of  Rotating  Stall 
Phenomena  in  Axial  Flow  Corrpressors , " with  Mr.  Marvin  A.  Stibich, 
AFAPL/TBC , as  Project  Engineer.  Dr.  Gary  R.  Ludwig  of  the  Cal span 
Corporation  was  technically  responsible  for  the  work.  Other  Calspan 
personnel  were:  Joseph  P.  Nenni,  John  C.  Erickson,  John  A.  Lordi, 

Gregory  F.  Homicz,  and  Rudy  H.  Arenat. 
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ABSTRACT 


This  report  presents  the  results  of  a research  program  that  had  two 
major  objectives.  The  first  objective  was  the  development  of  a prototype 
rotating  stall  control  system  which  was  tested  both  on  a low  speed  rig  and  a 
J-85-S  engine.  The  second  objective  was  to  perform  fundamental  studies  of  the 
flow  mechanisms  that  produce  rotating  stall,  surge  and  noise  in  axial  flow 
compressors  and  thereby  obtain  an  understanding  of  these  phenomena  that  would 
aid  attaining  the  first  objective.  The  work  is  reported  in  three  separate 
volumes.  Volume  I covers  the  fundamental  theoretical  and  experimental  studies 
of  rotating  stall;  Volume  II  covers  the  theoretical  and  experimental  studies  of 
discrete- tone  aerodynamic  noise  generation  mechanisms  in  axial  flow  compressors; 
and.  Volume  III  covers  the  development  and  testing  of  a prototype  rotating  stall 
control  system  on  both  the  low  speed  test  rig  and  the  J-85-S  engine. 

Volume  I describes  the  theoretical  and  experimental  investigation  of 
the  influence  of  distortion  on  the  inception  and  properties  of  rotating  stall 
for  an  isolated  rotor  row,  and  the  effects  of  close  coupling  of  a rotor  and 
stator  row  on  rotating  stall  inception.  The  experiments  were  conducted  in  the 
Calspan/Air  Force  Annular  Cascade  Facility,  which  is  a low  speed  compressor 
research  rig.  In  addition,  the  previously  developed  two  dimensional  stability 
theory  for  prediction  of  inception  conditions  was  extended  to  include  the 
effect  of  compressibility  and  the  development  of  a three  dimensional  theory 
was  in„  iated.  These  studies  led  to  the  following  key  results.  The  experi- 
mental studies  of  distortion  show  that  for  a single  blade  row  the  response  of 
the  blade  row  to  the  distortion  and  rotating  stall  are  uncoupled  phenomena  and 
may  be  explained  on  the  basis  of  a lineariced  analysis.  The  experimental 
studies  of  a closely  coupled  rotor-stator  pair  show  that  the  addition  of  a 
close’;'  spaced  stator  row  downstream  of  a rotor  row  delays  the  onset  of  rotating 
stall.  Moreover,  the  corresponding  theoretical  analysis  predicts  this  trend 
although  Muantitative  agreement  is  hampered  by  the  lack  of  appropriate  steady- 
state  loss  and  turning  performance  for  each  blade  row  at  the  required  operating 
conditions.  The  theoretical  investigation  of  the  effects  of  compressibility 
for  wholly  subsonic  flows  outside  the  blade  rows  indicates  that  the  effects  of 
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compressibility  do  not  alter  the  mechanisms  of  rotating  stall  as  deduced  from 
the  incomnressible  theory  in  that  the  rate  of  change  of  the  steady  state  loss 
curve  with  inlet  swirl  is  the  dominant  blade  row  characteristic  affecting  its 
stability.  Therefore,  if  the  steady  state  losses  are  known  for  the  compressible 
flow  condition,  the  linearized  stability  analysis  is  expected  to  apply. 

Volume  II  describes  a theoretical  and  experimental  study  of  discrete- 
tone  noise  generation  by  the  interaction  of  a rotor  and  a stator,  a.id  the  de- 
velopment of  a direct  lifting  surface  theory  for  an  isolated  rotor.  An  approxi- 
mate node!  has  been  developed  to  predict  the  sound  pressure  level  and  total 
power  radiated  at  harmonics  of  the  blade  passage  frequency  for  a rotor-stator 
stage.  The  analysis  matches  the  duct  acoustic  modes  for  an  annular  duct  with 
an  approximate  representation  of  the  unsteady  blade  forces  which  includes  com- 
pressibility effects.  Measurements  were  made  of  the  sound  pressure  levels 
produced  on  the  duct  wall  of  the  annular  cascade  facility  by  a rotor-stator 
pair.  Predictions  which  indicated  that  only  the  fourth  and  higher  harmonics 
could  be  excited  at  conditions  achievable  in  the  facility,  were  borne  out  by 
the  experiments.  The  calculations  of  the  sound  pressure  levels  for  the  propa- 
gating modes  were  significantly  below  the  measured  values.  This  discrepancy 
is  believed  to  result  from  inaccuracies  in  existing  models  of  rotor  wake 
velocity  profiles,  which  are  shown  to  have  a strong  influence  on  predictions 
of  the  sound  pressure  levels  of  the  higher  harmonics.  Volume  II  also  contains 
the  formulation  of  a direct  lifting  surface  theory  for  the  compressible, 
three-dimensional  flow  through  a rotor  row  in  an  infinitely  long  annular  duct. 

A detailed  derivation  is  given  for  the  linearized  equations  and  the  corres- 
ponding solutions  for  the  blade  t .lckness  and  loading  contributions  to  the 
rotor  flow  field.  The  governing  integral  equation  for  the  blade  loading  in  a 
lifting  surface  theory  is  obtained  for  subsonic  flow  and  progress  on  its 
solution  is  reported. 

Volume  III  describes  the  development  and  testing  of  3 prototype  ro- 
tating stall  control  system.  The  control  system  was  tested  on  the  low  speed 
compressor  research  rig  and  on  a J-85-5  turbojet  engine.  On  the  low  speed 
research  compressor,  the  control  was  tested  in  the  presence  of  circumferential 
inlet  distortion.  These  tests  were  performed  to  demonstrate  the  ability  of 
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tne  control  co  operate  satisfactorily  in  the  presence  of  inlet  distortion  and 
to  aid  in  the  selection  of  stall  sensor  configurations  for  the  subsequent 
engine  tests.  The  control  system  was  then  installed  on  a J-85-5  jet  engine  and 
its  performance  was  tested  under  sea  level  static  conditions,  both  with  and 
without  inlet  distortion.  On  the  engine,  the  stall  control  was  installed  to 
override  the  normal  operating  schedule  of  the  compressor  bleed  doors  and  inlet 
guide  vanes.  The  J-85-5  was  stalled  in  two  ways,  first  by  closing  the  bleed 
doors  at  constant  engine  speed,  and  second  by  decelerating  the  engine  with  the 
bleed  doors  partially  closed  at  the  beginning  of  the  deceleration.  \ total  of 
41  compressor  stalls  were  recorded  at  corrected  engine  speeds  between  48  and 
72  percent  of  the  rated  speed.  In  all  cases,  the  control  took  successful 
remedial  action  which  limited  the  duration  of  the  stall  to  525  milliseconds 
or  less. 
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SECTION  I 


INTRODUCTION 


The  useful  operating  range  of  a turbine  engine  compressor  is  greatly 
influenced  by  its  stalling  characteristics.  The  optimum  performance  of  a 
turbo-propulsion  system  is  usually  achieved  when  the  compressor  is  operating 
near  its  maximum  pressure  ratio  However,  this  optimum  is  generally  not 
attainable  because  it  occurs  close  to  compressor  stall  and  unstable  flow  con- 
ditions. Because  of  the  serious  mechanical  damage  that  may  result  during  com- 
pressor stall  cycles,  a factor  of  safety  (/stall  margin)  must  be  provided 
between  the  compressor  operating  line  and  the  stall  boundary.  This  is  usually 
done  by  preschvUul ing  the  primary  engine  controls.  However,  the  preschedu 1 i ng 
approach  can  lead  to  the  requirement  for  a large  stall  margin  in  order  to  keep 
the  engi  ie  from  stalling  under  all  possible  transient  and  steady  state  flight 
conditions.  It  is  clear,  then,  that  in  engine  control  system  that  can  sense 
incipient  destructive  unsteady  flow  in  a compressor  and  take  corrective  action 
would  allow  for  reduced  stall  margins  in  the  design  and  thus  lead  to  large 
engine  performance  and/or  efficiency  gains.  Recognition  of  this  fact  has  boon 
the  motivation  for  a continuing  program  of  research  that  the  AI-'.M’l.  has  sponsored 
at  Cal  span  dating  hack  to  HHiJ. 

The  work  at  Calspan  has  been  both  theoretical  and  experimental  in 
nature  and  has  been  aimed  at  obtaining  a sufficient  understand i ng  of  the  rot  at 
ing  stall  phenomena  such  that  its  onset  and  its  properties  can  he  predicted 
and  controlled.  The  capability  of  predicting  the  onset  of  rotating  stall  on 
isolated  blade  rows  of  high  huh  to  tip  ratios  in  low  speed  flows  was  demon 
strated  in  Reference  1.  In  addition,  the  basic  feasibility  of  developing  a 
rotating  stall  control  system  was  demonstrated  in  the  Calspan/Air  force  Annular 
Cascade  facility.  This  present  report  summarizes  the  latest  three  Near- 
research  program  at  Calspan.  The  specific  goals  of  the  present  program  were 
to  extend  the  fundamental  studies  of  rotating  stall  to  consider  the  effects  of 
compressibility,  blad _•  row  interaction  and  inlet  distortion;  and  to  extend  the 
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fundamental  aerodynamic  and  acoustic  analysis  of  flow  through  a compressor. 

In  addition,  the  rotating  stall  control  system  was  validated  by  successful 
ground  tests  on  a J-8S-5  turbojet  engine. 

The  work  is  reported  in  three  separate  volumes.  Volume  I entitled, 
"Basic  Studies  of  Rotating  Stall",  covers  the  theoretical  and  experimental 
work  on  the  effects  of  distortion  and  close  coupling  of  blade  rows  on  rotating 
stall  inceptio"  and  properties.  In  addition,  the  theoretical  analysis  of  com- 
pressibility is  treated  in  the  two-dimensior.a 1 approximation  and  the  initial 
development  of  a three-dimensional  theory  is  given.  Volume  II  entitled, 
"Investigation  of  Rotor-Stator  Interaction  Noise  and  Lifting  Surface  Theory 
for  a Rotor",  describes  the  development  of  a lineariced  lifting  surface  theory 
for  the  subsonic  compressible  flow  through  an  isolated  rotor  row.  In  addition, 
a theoretical  and  experimental  study  of  the  noise  generated  by  the  interaction 
of  a rotor  and  stator  is  described.  Volume  III  entitled,  "Development  of  a 
Rotating  Stall  Control  System",  describes  the  development  and  testing  of  the 
control  system  installed  on  a low  speed  research  compressor  and  on  a J-S5-S 
turbojet  engine. 

Volume  I has  been  divided  into  two  main  sections  which  are  Experimental 
Rotating  Stall  Research  and  Theoretical  Rotating  Stall  Research.  The  section 
on  experimental  research  describes  the  low  speed  tests  run  in  the  Calspan/Air 
Force  Annular  Cascade  Facility  to  investigate  the  effects  of  inlet  distortion 
and  blade  row  interference  on  routing  stall.  The  section  on  theoretical  re- 
search describes  the  two  dimensional  stability  theory  that  was  developed  for 
correlation  with  the  above  experimental  data.  In  addition  this  section  con- 
tains the  summary  of  the  results  of  extension  of  the  two  dimensional  theory  to 
compressible  flow  and  the  initial  developments  of  a three  dimensional  theory. 
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bcCTION  II 

EXPERIMENTAL  ROTATING  STALL  RESEARCH 


As  a part  of  the  work  under  a previous  program  Contract  AF  33(615) -55S7 
an  annular  cascade  facility  was  designed  and  fabricated.  Its  principal  purpose 
is  to  provide  detailed  fundamental  experimental  data  during  and  prior  to  the 
occurrence  of  rotating  stall  in  order  to  improve  our  understanding  of  the  phe- 
nomena and  for  use  as  a guide  in  improving  the  theoretical  analysis.  The 
facility  has  also  been  used  to  evaluate  the  operation  of  a prototype  rotating 
stall  control  system  and  to  provide  acoustic  data  for  comparison  with  theory. 

The  acoustic  measurements  are  described  in  Volume  II  of  this  report  and  the 
control  system  tests  are  described  in  Volume  III.  This  section  presents  the 
fundamental  experiments  on  rotating  stall. 

In  previous  programs  using  the  annular  cascade  facility,  rotating 
stall  lias  been  induced  on  a variety  of  stator  rows  with  differing  geometries 
and  on  a rotor  row  at  various  stagger  angles.  The  experimental  inception  con- 
ditions and  properties  of  rotating  stall  were  measured.  All  of  these  investi- 
gations were  conducted  in  the  absence  of  deliberate  inlet  distortion.  We  at- 
tempted to  minimize  distortion  although  some  unavoidable  radial  distortion  was 
present  because  of  the  use  of  inlet  guide  vanes  in  most  test  configurations. 

The  distortion  was  measured  but  no  attempt  was  made  to  vary  or  investigate  the 
distortion.  Since  it  is  known  that  inlet  distortion  can  influence  surge  and 
rotating  stall  inception,  a major  portion  of  the  current  program  was  devoted 
to  investigating  the  effects  of  circumferential  inlet  distortion  on  rotating 
stall.  This  study  consisted  of  two  parts.  In  the  first  part,  an  isolated 
rotor  was  tested  in  the  presence  of  distortion  patterns  which  wen;  stationary 
in  the  laboratory  frame  of  reference.  In  the  second  part,  one  of  the  distortion 
patterns  was  rotated  about  the  compressor  axis  at  a series  of  angular  velocities 
while  the  rotor  velocity  was  held  fixed.  The  latter  studv  was  designed  to 
determine  if  the  time  of  passage  of  a rot  ir  blade  through  a distortion  cell 
plays  a significant  role  in  blade  row  sensitivity  to  inlet  distortion. 

In  addition  to  the  inlet  distortion  studies,  the  current  program  in- 
cluded a limited  examination  of  the  effects  of  rotor-stator  interference  on 
the  inc  ption  of  rotating  stall.  The  rotor  row  and  stator  row  chosen  for  these 
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experiments  had  previously  been  tested  in  isolation  so  that  the  conditions 
under  which  the  individual  rotor  row  and  stator  row  undergo  rotating  stall 
inception  were  known.  Thus  the  presence  of  interference  effects  between  the 
two  rows  when  combined  as  a rotor-stator  stage  could  be  determined  by  compari- 
son with  the  individual  blade  row  performances.  The  results  were  also  used 
for  comparison  with  the  predictions  of  rotating  stall  theory  for  multiple 
blade  rows. 

The  order  of  presentation  of  the  experimental  program  in  the  remainder 
of  this  section  is  as  follows.  A description  of  the  annular  cascade  facility 
is  presented  in  Section  1 1 - A . Section  1 1 - B presents  the  results  of  the  study 
of  rotor-stator  interference  effects  on  rotating  stall.  Section  II-C  presents 
the  results  of  the  rotating  stall  experiments  in  the  presence  of  inlet  dis- 
tortion. This  section  includes  the  work  on  both  stationary  inlet  distortion 
and  on  rotating  inlet  distortion.  Finally,  in  Section  II-D,  a summary  of  the 
experimental  program  un  rotating  stall  is  presented  along  with  those  vonclu- 
siens  which  can  be  drawn  from  inspection  of  the  results. 

A.  DESCRIPTION  OF  ANNULAR  CASCADE  FACILITY 

A detailed  description  of  the  arnular  cascade  facility  has  been  pre- 
sented in  Reference  1.  Rather  than  repeat  much  of  this  detail,  the  following 
description  is  oriented  towards  noting  changes  in  the  facility  nvde  during 
the  current  program. 

1 . General  Description 

The  annular  cascade  facility  consists  of  a test  section  built  around 
the  outer  front  casing  of  a J-79  jet-engine  compressor  with  a Calspan  fabri- 
cated hub.  The  facility  includes  a bell-mouth  inle*  on  the  outer  casing  and 
a bullet  nose  on  the  hub  to  provide  a smooth  flow  of  air  to  the  test  section. 
Outlet  ducting  is  connected  to  an  independently  variable  source  of  suction  to 
provide  the  required  flow  through  the  annulus.  An  electrically  powered  two- 
speed  axial  flow  fan  is  used  as  the  source  of  suction.  Continuous  control  of 
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the  mass  flow  is  achieved  through  the  use  of  variable  inlet  guide  vanes  to 
the  fan  and  a variable  damper  in  the  fan  exit  flow.  Photographs  of  the  com- 
plete annular  cascade  facility  and  of  the  test-section  portion  of  the  facility 
are  shown  in  Figures  1 and  2,  respectively. 

The  test  section  of  the  annular  cascade  forms  a circular  annulus 
with  an  outer  diameter  of  29. 35  inches  and  an  inner  diameter  of  23.35  inches 
which  provides  a hub-to-tip  ratio  of  0.80.  The  outer  casing  will  accept  up 
to  six  variable  stagger  angle  stator  rows.  The  hub  has  provision  for  two 
rotor  rows  at  the  third-  and  fifth-stage  rotor  locations  of  the  J-79  compressor. 
At  present,  speed  control  on  each  rotor  hub  is  independently  variable  in 
either  direction  of  rotation.  The  overall  combination  of  variable  stagger 
angle  stators,  independently  powered  rotors,  and  independent  mass  flow  control 
provides  exceptional  versatility  to  the  complete  test  rig. 

2 . Rotating  Hub  Drive  Systems 

During  the  course  of  the  current  investigations,  two  different  drive 
systems  were  used  to  power  the  rotating  hubs  in  the  annular  cascade  facility. 

The  original  drive  system  allowed  both  rotors  to  be  driven  at  the  same  speed 

or  one  rotor  to  be  driven  while  the  other  was  held  fixed.  All  experiments 

which  required  the  use  of  only  one  powered  rotor  were  performed  with  this 
drive  system.  However,  the  experiments  with  rotating  distortion  required  each 
rotor  hub  to  have  independent  speed  control.  Thus  the  rotor  drive  system  was 
modified  to  provide  this  capability.  The  drive  system  modification  did  not 
alter  any  dimensions  within  the  circular  annulus  test  section  of  the  facility. 
Both  rotor  drive  systems  are  described  in  the  following  paragraphs. 

Details  of  the  original  rotating  hub  installation  in  the  annular 
cascade  are  shown  in  Figure  3.  The  main  features  of  the  design  are  as  follows. 
There  are  two  rotor  assemblies,  one  centered  at  the  axial  location  of  the 
third  stage  rotor  in  the  J-79  compressor  and  one  centered  at  the  J-79  fifth 
stage  rotor.  Both  of  these  rotors  are  driven  by  a common  drive  shaft.  How- 
ever, either  rotor  can  be  decoupled  from  the  drive  shaft  and  held  stationary 
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while  the  other  rotor  is  driven.  Alternatively,  both  rotors  can  be  driven, 
but  only  at  identical  velocities  in  the  same  direction.  Only  the  rotor  as- 
semblies rotate,  the  outer  skin  on  the  hub  upstream  and  downstream  of  the 
rotors  and  between  the  rotors  is  held  stationary.  The  drive  shaft  to  the 
rotors  is  powered  by  a 24  horsepower  hydraulic  motor.  Rotational  speed  is 
infinitely  variable  in  either  direction  between  zero  and  approximately  1500  rpm. 
An  external  hydraulic  pump  system  powered  by  a 'V'1  horsepower  electric  motor 
is  used  to  provide  power  for  the  hydraulic  motor 

The  modified  rotating  hub  installation  is  shown  m Figure  4.  In 
this  design  both  the  third  and  fifth  stage  rotor  hubs  remain  in  the  same 
axial  locations  as  in  Figure  3.  The  third  stage  : otor  shaft  is  driven  directly 
by  the  original  drive  shaft,  hydraulic  motor,  and  hydraulic  pump  system. 

The  fifth  stage  rotor  is  mounted  on  a concentric  drive  shaft  which  is  chain 
driven  by  a second  hydraulic  motor.  Power  for  this  second  hydraulic  motor  is 
provided  by  ?.  separate  and  larger  hydraulic  pump  system  with  a 50  horsepower 
electric  motor.  With  this  modified  hub,  both  rotor  hubs  have  completely  inde- 
pendent speed  controls  and  can  be  rotated  in  either  direction.  Here  again, 
only  the  rotor  assemblies  rotate,  the  outer  skin  on  the  hub  upstream  and 
downstream  of  the  rotors  and  between  the  rotors  is  held  stationary.  Maximum 
speed  capability  on  the  third  stage  rotor  hub  is  approximately  1500  rpm.  Max- 
imum capability  on  the  fifth  stage  rotor  hub  has  nut  been  tested,  but 

it  is  _,er  than  that  of  the  third  stage  rotor  because  of  the  larger  hydraulic 
pump  system. 

The  two  rotor  hubs  have  provision  for  accepting  different  numbers  of 
ro*;or  blades.  The  stage  3 rotor  hub  i-  designed  to  accept  3b  blades  while 
the  stage  5 rotor  hub  is  designed  to  accept  16  blades.  The  method  of  fastening 
the  blades  to  either  hub  allows  individual  adjustment  of  the  stagger  angle. 
Details  of  the  blade  fastenings  are  shown  m an  ir<<  rt  in  Figure 


In  addition  to  rotor  Males , two  symmetric  multi -lobe 
creetis  were  designed  for  mounting  on  the  third,  stags  ro*or  h ib . 
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rotated  about  the  compressor  axis.  The  screens  are  perforated  steel  plates. 

The  mounting  system  uses  the  same  holes  as  the  rotor  blades.  Details  of  the 
screen  mounting  system  are  shown  in  an  insert  in  Figure  4.  This  system 
allows  configuration  changes  to  be  made  without  disassembly  of  the  hub  portion 
of  the  rotating  rig. 

3 . Rotor  and  Stator  Blading 

The  tests  performed  during  this  program  involved  either  an  isolated 
rotor  in  the  presence  of  inlet  distortion  or  a closely  coupled  rotor-stator 
stage.  The  fifth  stage  rotor  and  stator  blades  from  a J-"9  compressor  were 
used  in  these  tests,  unmodified  except  for  Made  length  which  was  shortened 
in  order  to  fit  ir.to  the  three-inch  annulus  of  the  test  rig.  In  each  case, 
the  spanwise  portion  of  the  blades  closest  to  the  outer  compressor  casing  was 
used.  The  set  of  rotor  blades  has  been  designated  as  Rotor  Set  No.  1 and  the 
set  of  stator  blades  as  Stater  Set  No.  ! in  keening  with  the  nomenclature  used 
in  Reference  1.  The  geometric  character ist ics  ot  these  blades  have  been 
presented  in  Reference  1 along  with  detailed,  measurements  ot  tun  ing  aim  loss 
performance  when  used  as  single  blade  rows  m isolation.  The  extensile  per- 
formance data  are  not  repeated  herein  except  where  the;.’  are  used  for  purposes 
of  comparison  wi*h  the  current  test  results.  The  blade  geometric  character- 
istic' are  repeated  below. 

Rotor  'set  No.  1 has  the  sine  -umber  <f  blades  and  i«  at  the  same 
axial  location  as  -he  fifth,  stage  rot  ir  in.  *ue  enigma’.  l-'D  compressor.  These 
blades  have  a N.V*A  *>*  *•»  *h  • cr*  • «s  '■  .*  > > nut  ion  on  a circular  a rc  m-ai 

line. 
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TABLE  I 


Geometric 

Blade  Length 
Blade  Chord 
Blade  Thickness; 


Characteristics 


Outer  Diameter 
Mid-Annulus 
Inner  Diameter 


Rotor  Set  N’o.  1 

2.98 

1.448 

0.072 

0.086 

0.113 


inches 

inches  (constant) 

inch 

inch 

inch 


Blade  Camber  Angle  (angle  between  tangents 
to  mean  camber  line  of  leading  and  trailing 


edges)  ; 

Outer  Diameter 

23.5  degrees 

Mid-Annulus 

28.0  degrees 

Inner  Diameter 

32.7  degrees 

Blade  Twist 

See  Below 

Number  of  Blades 

46 

Solidity  at  Mid- 

Annulus 

0.805 

S 


The  blades  in  Rotor  Set  No.  1.  have  a nominal  twist  of  14. 1 degrees 
over  tiie  three  inch  span.  However,  measurements  showed  that  in  practice,  the 
twist  over  the  center  2-4  inches  in  these  production  J-79  blades  varied  be- 
tween 10.8  and  14.4  degrees  with  an  average  value  of  12.9  degrees.  The  blades 
were  installed  in  a sequence  which  distributed  the  nonuniformities  in  twist 
over  the  complete  circumference. 

A blade  tip  clearance  problem  was  encountered  on  installation  of 
the  rotor  assembly  in  the  annular  cascade.  The  J-79  compressor  outer  casing 
used  in  the  annular  cascade  is  about  one-tenth  inch  out  of  round  with  the 
split  line  diameter  being  the  largest.  In  final  assembly,  the  rotor  was 
mounted  in  the  casing  with  shims  along  the  split  line  which  provided  the 
following  tip  clearances:  0.030  inch  on  the  top  and  bottom  and  0.04n  to  0.048 

inch  on  the  sides.  Since  the  stagger  angles  of  the  ro^or  blades  are  adjustable 
the  quoted  tip  clearances  apply  only  to  the  mid-chord  pivot  points.  Clearance 
at  the  blade  leading  and  trailing  edges  vary  slightly  with  changes  in  stagger 
angle  from  the  angle  for  which  the  tip  contour  was  determined.  The  stagger 
angle  used  to  determine  the  tip  contour  was  the  maximum  which  was  expected  to 
be  used.  The  choice  of  lower  reference  stagger  angles  would  have  resulted  in 
interference  between  the  blades  and  outer  casing  at  stagger  angles  much  larger 
than  the  reference. 

Stator  Set  No.  1 has  the  same  number  of  blades  and  is  located  at  the 
same  axial  location  as  the  fifth  stage  stator  row  in  the  original  J-"9  com- 
pre.'.sor.  Its  geometric  characteristics  are  listed  in  I \B!.i  ii. 

1 . Instrumentation 


The  outer  casing  of  the  annular  cascade  test  section  has  been 
modified  to  allow  circumferential  and  radial  traverses  with  hot-wire  and 
total  pressure  probes  upstream  and  downstream  of  the  fifth  stage  location. 

These  traverse  locations  are  shown  in  Figure  2.  Two  probes  can  be  accommodated 
simultaneously  in  either  of  these  traverse  mounts  with  adjustable  i n nr: 
ferential  spacing  between  the  two  probes.  In  addition,  simple  radial  traverse- 
can  be  made  at  various  other  locations  in  the  test  section. 


TABLE  II 


Geometric  Characteristics  of  Stator  Set  No. 


Blade  Length 

Blade  Chord;  Outer  Diameter 
Inner  Diameter 

Blade  Thickness;  Outer  Diameter 
Inner  Diameter 

Blade  Camber  Angle  (angle  between  tangents 
to  mean  camber  line  of  leading  and  trailing 
edges') 


Outer  Diameter 
Blade  Twist;  Mid-Annulus 
positive  twist  Inner  Diameter 
reduces  stagger 

Number  of  Blades 

Solidity  at  Mid-Annulus 


in 


1 

2.98  inches 

1.316  inches 
1.290  inches 

0.15<.’  inch 
0.1132  inch 

35.6  degrees 

0 degree 
0.28  degree 
1.42  degree 

54 

0.85 


J 


A linearized  two-channel,  constant -temperature , hot-wire  anemometer 
system  was  used  in  conjunction  with  a crossed-wire  probe  fcr  the  velocity  and 
swirl  angle  measurements.  Readout  for  the  hot-wire  system  was  made  by  means 
of  a two-channel  integrator-digital  voltmeter  system.  Each  linearized  hot- 
wire signal  was  integrated  for  60  seconds  and  the  time  average  was  calculated 
from  the  readings  on  the  digital  voltmeters. 

Total  pressure  surveys  were  made  through  circumferential  traverses 
with  a multiple-tube  total  pressure  rake.  The  total  pressure  rake  was  aligned 
with  the  flow  by  using  a Conrad  arrowhead  style  yawmeter  incorporated  on  the 
rake.  In  all  tests,  the  pressures  detected  by  the  rake  were  photographically 
recorded  from  a multitube  inclined  manometer. 

In  most  cases,  rotating  stall  was  detected  and  the  number  of  cells 
and  propagation  velocity  were  measured  by  using  outer-wall  static  pressure 
taps  connected  to  pressure  transducers.  Two  pressure  taps  at  different  cir- 
cumferential locations  and  at  an  axial  location  corresponding  to  the  quarter- 
chord  on  the  rotor  blades  were  used.  In  the  tests  with  rotating  inlet  dis- 
tortion, rotating  stall  was  detected  in  several  ways.  The  most  informative 
method  consisted  of  measuring  the  mid-annulus  total  pressure  fluctuations 
just  upstream  of  the  rotor. 

In  all  of  the  experiments,  mass  flow  through  the  test  section  was 
determined  by  measuring  the  dynamic  pressure  in  the  constant -area  annulus 
upstream  of  the  test  section. 


B.  INVESTIGATION  OF  EFFECT  OF  ROTOR- STATOR  INTERFERENCE  ON- 

ROTATING  STALL 

Previous  studies  of  rotating  stall  in  the  annular  cascade  facility 
were  concentrated  on  isolated  rows  of  rotors  on  stators.  Although  the  stator 
row  studies  required  the  use  of  upstream  guide  vanes  to  generate  the  swirlea 
flow  necessary  to  provide  the  conditions  for  rotating  stall,  the  axial 
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separation  between  the  guide  vane  row  and  stator  row  was  large  compared  to 
the  typical  spacing  of  a rotor-stator  pair  usually  found  in  a compressor 
stage.  Moreover,  in  most  of  these  tests,  both  the  guide  vane  row  and  the 
stator  row  were  stationary.  Thus  the  current  experiments  were  undertaken  to 
determine  if  interference  between  a closely  coupled  rotor-stator  pair  changes 
the  conditions  for  inception  of  rotating  stall  from  those  which  are  found  on 
isolated  blade  rows. 

The  notation  for  these  experiments  is  shown  in  Figure  5 and  a sketch 
of  the  annular  cascade  with  its  geometric  dimensions  is  shown  in  Figure  6. 

The  rotor-stator  stage  used  in  the  annular  cascade  consisted  of  Rotor  Set  No.  1 
followed  by  Stator  Set  No.  1.  These  are  modified  blade  rows  from  the  fifth 
stage  of  a J-79  compressor  (See  Section  II. A. 3).  The  steady  state  performance 
of  these  blade  rows  had  been  determined  previously  (Reference  1)  for  the  con- 
ditions near  rotating  stall  inception  on  each  of  these  blade  rows  in  isolation. 
In  addition,  the  inception  point  and  properties  of  rotating  stall  of  each 
blade  row  in  isolation  had  also  been  determined.  In  the  current  work,  the 
inception  point  and  properties  of  rotating  stall  on  the  combined  blade  rows 
were  studied.  The  techniques  used  to  measure  rotating  stall  properties  are 
the  same  as  those  used  on  the  isolated  rotor.  These  have  been  described  in 
Reference  1. 


The  tests  were  performed  for  one  stagger  angle  setting  of  the  rotor 
blades  ( <5»m  = 40  degrees)  and  for  four  stagger  angle  settings  of  the  stator 
blades  ( 6 iH  = 28.2,  37.2,  47.2  and  57.2  degrees).  Since  both  rotor  and  stator 
blades  are  twisted,  the  quoted  stagger  angles,  6nM  and  <JiH  , correspond  to 
the  values  at  mid-annulus  in  the  onnular  cascade.  In  addition  to  determining 
the  rotating  stall  inception  point,  propagation  velocity,  and  number  of  cells, 
the  static  pressure  rise  across  the  stage  was  also  measured  over  a range  of 
rotor  rpm's  for  each  of  the  stator  stagger  angle  settings. 

A summary  of  the  inception  points  found  in  these  tests  is  presented 
in  Table  III  along  with  the  isolated  rotor  results  for  comparison.  In  all 
cases  tested,  the  presence  of  rotating  stall  was  detectable  from  signals 


NOTE:  MEAN  AXIAL  VELOCITY,  U = 60  fps  FOR  ALL  TESTS. 


generated  by  pressure  transducers  installed  near  the  blade  tips  of  both  the 
rotor  row  and  the  stator  row.  Thus,  it  is  not  possible  to  determine  absolutely 
which  blade  row  was  responsible  for  rotating  stall  when  it  occurred.  However, 
for  the  two  lowest  stator  stagger  angles,  diM  a 28.2  and  37.2  degrees,  and  to 
a lesser  extent  for  6tH  = 47.2  degrees,  the  stall  signals  displayed  character- 
istics similar  to  rotating  stall  on  the  isolated  rotor.  At  the  highest  stator 
stagger  angle,  the  rotating  stall  signals  were  very  irregular  on  both  blade 
rows  but  the  general  features  looked  more  like  those  associated  with  stator 
stall.  In  any  event,  both  blade  rows  in  the  closely  coupled  stage  appeared 
to  participate  strongly  in  the  rotating  stall  phenomenon. 

In  Table  III,  the  rotor  speed  at  inception  with  increasing  rotor 
speed  is  shown  along  with  the  corresponding  overall  inlet  swirl  angle  /§Q|t 
relative  to  the  rotor.  Similar  data  are  shown  for  the  point  at  w.  ich  rotating 
stall  dies  out  as  rotor  speed  i-  decreased.  The  sign  convention  for  angles  is 
given  in  Figure  5.  For  all  of  the  cases  tested,  the  addition  of  a stator 
row  downstream  of  the  rotor  row  delayed  the  inception  of  rotating  stall.  The 
delay  in  inception  increased  as  the  stator  stagger  angle  was  increased.  In 
terms  of  rotor  speed,  the  delay  in  inception  was  substantial  (»25^)  at  the 
two  highest  stator  stagger  angles.  The  case  with  dSM  = 57.2  degrees  is  probablv 
not  of  practical  interest  since  the  stator  blades  were  either  very  lightly 
loaded  or  even  negatively  loaded.  Nevertheless,  these  data  could  be  useful 
for  assessing  the  validity  of  the  two  blade  row  rotating  stall  theory  under 
widely  varying  conditions. 

The  delay  in  rotating  stall  inception  to  higher  rotor  speeds  on  the 
coupled  rotor-stator  stage  has  a large  effect  on  the  static  pressure  rise  at- 
tainable by  the  stage.  Moreover,  the  behavior  of  the  pressure  rise  near  stall 
appears  to  be  dependent  on  the  degree  of  loading  of  the  stator  blades.  These 
effects  are  discussed  in  the  following  paragraphs. 

The  static  pressure  rise  data  measured  across  the  isolated  rotor 
and  across  the  rotor-stator  stage  are  shown  in  Figures  7 and  8 respectively. 

The  data  were  measured  both  at  the  hub  and  at  the  tip  through  pressure  taps 
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located  on  the  inner  and  outer  casings  of  the  annular  cascade.  In  these  ex- 
periments, the  pressure  rise  data  were  obtained  over  a range  of  rotor  speeds 
with  the  mean  axial  flow  velocity,  U,  , held  constant  at  60  fps.  The  pressure 
rise  has  been  made  dimensionless  by  dividing  by  the  dynamic  pressure  correspond- 
ing to  this  value  of  U,  . On  each  figure  the  presence  of  rotating  stall  is 
indicated  by  tails  on  the  symbols  and  the  rotor  speed  at  rotating  stall  in- 
ception is  indicated.  In  instances  where  inception  displayed  hysteresis, 
data  are  shown  both  with  and  without  the  presence  of  rotating  stall. 

As  expected,  the  isolated  rotor  data  in  Figure  7 show  that  the  static 
pressure  rise  across  the  blade  row  is  greater  at  the  tip  than  at  the  hub,  and 
that  the  difference  increases  with  increasing  rotor  speed  in  the  unstalled 
region.  This  is  in  contrast  with  the  pressure  rise  characteristic  for  a 
stator  row  which  generally  displays  the  greatest  static  pressure  rise  at  the 
hub.  See  for  example  Figure  28  of  Reference  2.  When  the  rotor  and  stator  ■’re 
combined  to  make  a stage,  as  in  Figures  8(a)  through  8(d),  the  difference 
between  hub  and  tip  static  pressure  rise  on  the  rotor  and  stator  tend  to 
cancel  provided  that  the  two  blade  rows  are  well  matched.  In  these  figures, 
the  data  iur  diM  * 28.2  degrees  provided  the  closest  match  between  rotor  and 
stator  characteristics.  As  the  stator  stagger  angle  is  increased,  the  blade 
row  matching  deteriorates  and  the  difference  between  the  hub  and  tip  static 
pressure  rise  increases.  The  data  in  Figure  8(d)  for  dSM  = S7.2  degrees  is 
a particularly  bad  mismatch.  The  data  for  the  isolated  rotor  displays  a 
greater  pressure  rise  over  the  full  range  of  rotor  speeds  tested  than  the 
data  for  the  stage  with  <SiH  = 57.2  degrees.  This  is  true  even  for  the  rotor 
speed  range  whe *e  the  isolated  rotor  was  in  rotating  stall  while  the  stage 
was  not. 


The  stage  pressure  rise  data  behave  as  expected  in  the  rotor  speed 
range  where  rotating  stall  is  absent.  At  any  given  rot^r  speed,  the  highest 
pressure  rise  is  attained  when  the  stator  row  is  most  highly  loaded 
( rfSM  = 28.2  degrees)  and  the  pressure  rise  decreases  progressively  as  the 
stator  row  is  unloaded  ( (J>M  increasing).  However,  the  pressure  rise  data 
show  an  interesting  feature  as  inception  is  approached.  The  attainable  pres- 
sure rise  just  before  inception  increases  in  value  as  the  stator  row  is  un- 
loaded (except  for  the  badly  mismatched  ruse  = 57.2  degrees).  Apparently 
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decreasing  the  loading  on  the  stator  row  increases  the  rotor  speed  at  incep- 
tion by  an  amount  large  enough  that  the  maximum  attainable  pressure  rise  is 
increased  before  rotating  stall  occurs. 

With  the  isolated  rotor  (Figure  7)  a small  hysteresis  in  rotating 
stall  inception  was  observed  and  the  pressure  rise  curves  display  a small 
discontinuity  in  both  magnitude  and  slope  in  the  hysteretic  region.  For  the 
stage  with  d,M  ■ 28.2  and  37.2  degrees  (Figures  8(a)  and  8(b))  rotating  stall 
inception  was  again  hysteretic  and  the  hysteresis  was  not  large  in  terms  of 
rotor  speed.  However,  the  discontinuity  in  the  pressure  rise  curves  is  much 
more  pronounced  especially  for  6Sh)  * 28.2  degrees.  At  the  two  highest  stator 
stagger  angles  (Figures  8(c)  and  8(d))  no  hysteresis  was  observed  at  inception. 
In  these  cases  the  pressure  rise  curves  show  a discontinuity  in  slope  at  in- 
ception but  there  does  not  appear  to  be  a significant  discontinuity  in  mag- 
nitude. 


The  experimentally  determined  rotating  stall  propagation  velocities 
and  numbers  of  cells  are  shown  in  Figure  9(a)  for  stator  stagger  angles  of 
28.2  and  37.2  degrees  and  in  Figure  9(b)  for  stator  stagger  angles  of  47.2 
and  57.2  degrees.  Similar  data  measured  previously  on  the  isolated  rotor  are 
shown  on  each  figure  for  comparison.  At  the  three  lowest  stator  stagger 
angles  the  dimensionless  propagation  velocities  relative  to  the  rotor  row  are 
substantially  higher  than  these  found  on  the  isolated  rotor.  There  also  ap- 
pears to  be  a small  decrease  in  relative  propagation  velocity  as  the  stator 
stagger  angle  is  increased.  Near  inception,  a single  stall  cell  was  found  to 
be  propagating  for  stator  stagger  angles  of  28.2,  37.2  and  47.2  degrees  as 
well  as  for  the  isolated  rotor.  After  inception,  the  number  of  cells  remained 
at  one  for  d'SH  = 28.2  and  37.2  degrees  but  increased  to  a maximum  of  four  for 
ds H=  47.2  degrees.  The  number  of  stall  cells  on  the  isolated  rotor  increased 
to  two  well  after  inception.  Thus  for  these  cases,  the  number  of  cells  near 
inception  is  similar  to  the  isolated  rotor  and  after  inception  there  are  some 
small  differences.  The  characteristics  of  the  pressure  signals  from  the  ro- 
tating stall  sensors  looked  like  those  associated  with  rotating  stall  on  the 
isolated  rotor. 


At  a st'tor  stagger  angle  of  57.2  degrees,  the  rotating  stall 
properties  changed  considerably  from  those  found  with  the  lower  stator  stagger 
angles.  The  relative  propagation  velocities  dropped  to  values  near  those 
found  on  the  isolated  rotor  and  the  number  of  cells  increased  to  between  seven 
and  eight.  Moreover,  as  noted  previously,  the  pressure  signals  from  the  ro- 
tating stall  sensors  had  general  features  which  looked  more  like  those  asso- 
ciated with  stator  stall  than  with  rotor  stall. 

In  summary,  the  comparison  of  rotating  stall  properties  observed  on 
the  closely  coupled  stage  with  those  observed  on  the  isolated  rotor  indicates 
that  the  interaction  between  the  blade  rows  in  a stage  can  have  a substantial 
effect  on  the  inception  and  properties  of  rotating  stall.  The  most  interesting 
result  was  obtained  on  the  stage  with  the  stator  stagger  angle  set  at  47.2 
degrees.  With  this  configuration,  the  static  pressure  rise  measured  at  in- 
ception eliminated  the  sharp  drop  in  pressure  rise  when  rotating  stall  first 
occurred.  In  Section  III,  the  results  of  preliminary  correlations  between 
the  predictions  of  the  two-blade  row  theory  and  the  results  of  these  experi- 
ments are  presented.  These  correlations  are  encouraging  in  that  the  theory 
predicts  that  blade  row  coupling  does  stabilize  the  stage  and  appears  to  pro- 
vide an  accurate  quantitative  value  for  the  delay  in  inception  for  a stator 
stagger  angle  of  28.2  d<grees.  However,  the  theoretical  predictions  are  based 
on  large  extrapolations  of  the  stator  steady  state  turning  and  loss  performance. 
These  extrapolations  were  necessary  because  it  was  found  that  the  data  avail- 
able for  Stator  Set  No.  1 did  not  fall  within  the  range  of  inlet  conditions 
which  occur  when  the  stators  are  coupled  to  Rotor  Set  No.  1.  Moreover,  at 
the  larger  stator  stagger  angles  tested,  the  large  delay  in  inception  leads 
to  the  problem  of  estimating  the  steady  state  losses  on  the  rotor  under  con- 
dition-’ which  normally  induce  rotating  stall  on  the  isolated  rotor.  Thus  it 
is  planned  to  extend  the  steady  state  loss  and  turning  performance  measure 
ments  cn  the  isolated  stator  to  include  the  required  inlet  conditions  and  also 
to  measure  the  performance  of  the  rotor  row  and  stator  row  when  combined  in 
a closely  coupled  stage. 
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c. 


INVESTIGATION  OF  EFFECT  OF  STATIONARY  AND  ROTATING  INLET 
DISTORTION  ON  ROTATING  STALL 


There  have  been  numerous  tests  which  used  circumferential  distortion 
at  the  inlet  to  multistage  axial  flow  compressors.  These  tests  are  usually 
performed  in  order  to  determine  the  reduction  in  stall  margin  of  the  compres- 
sor when  subjected  to  such  distortion.  In  some  cases,  the  compressor  appears 
to  be  most  sensitive  to  a particular  circumferential  extent  of  the  distortion 
pattern.  In  fact,  there  is  evidence  that  multi-lobe  distortion  patterns,  all 
of  which  have  the  same  total  extent,  can  have  different  effects  on  compressor 
stall  margin  (e.g.,  References  3 and  4).  This  indicates  that  the  observed 
reductions  in  stall  margin  could  depend  on  a dii,  .-nsionless  reduced  frequency 
or  Strouhal  Number  which  is  formed  from  the  time  of  passage  of  a rotor  blade 
through  a distortion  pattern,  a typical  relative  flow  velocity,  and  a typical 
rotor  dimension  in  the  direction  of  the  relative  flow.  A stall-delay  concept 
based  on  the  above  idea  has  been  recommended  recently  by  Korn'*  for  the  im- 
provement of  parallel  compressor  theories. 

Investigations  of  the  above  type  are  very  useful  for  estimating  the 
effect  of  inlet  distortion  on  compressors  of  similar  design  and  for  providing 
some  understanding  of  the  parameters  involved.  However,  they  fall  short  of 
providing  an  understanding  which  is  sufficient  to  predict  the  way  in  whicli 
inlet  distortion  can  initiate  rotating  stall  and  surge.  For  example,  it  is 
possible  that  the  geometry  itself,  such  as  the  number  of  lobes,  can  effect 
rotating  stall  inception  independently  of  the  passage  time  of  a rotor  blade 
through  a distortion  lobe.  Geometry  dependent  distortion  pattern  effects 
similar  to  those  obtained  with  the  stall  delay  concept  have  been  proposed  in 
Reference  6 by  considering  only  the  effect  of  a diffuser  downstream  of  the 
compressor. 

The  work  presented  in  the  following  paragraphs  was  designed  to 
study  the  effect  of  circumferential  inlet  distortion  on  a simple  configuration 
consisting  of  an  isolated  rotor  in  a constant  area  annulus.  The  configuration 
was  kept  simple  in  order  to  facilitate  comparison  with  the  theoretical  studies 
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of  Section  III.  In  addition,  these  experiments  were  designed  to  search  for 
unambiguous  evidence  that  rotating  stall  inception  may  be  sensitive  to  a re- 
duced frequency  which  contains  the  time  of  passage  of  the  rotor  blades  through 
the  distortion  pattern.  In  order  to  investigate  reduced  frequency  effects, 
the  experimental  program  consisted  of  two  parts.  In  the  first  part,  the  ef- 
fect of  stationary  multi-lobe  distortion  patterns  on  rotating  stall  was  studied. 
Two  symmetric  multi- lobe  distortion  patterns  were  used.  Each  pattern  covered 
the  same  total  area  but  one  consisted  of  two  lobes  and  one  consisted  of  four 
lobes.  In  he  second  part,  the  two-lobe  dirtorticn  pattern  was  rotated  about 
the  compressor  axis  at  different  speeds  while  holding  the  rotor  speed  constant. 
The  latter  study  allowed  variations  in  the  relative  reduced  frequency  seen  by 
the  rotor  without  substantial  changes  in  the  geometry  of  the  distortion  pat- 
tern. Interpretation  of  changes  in  rotating  stall  inception  observed  in 
these  experiments  is  subject  to  less  ambiguity  than  evidence  inferred  indi- 
rectly from  multi- lobe  distortion  investigations  in  which  the  geometry  of  the 
distortion  pattern  is  changed. 

Both  the  stationary  and  rotating  distortion  investigations  were 
performed  with  Rotor  Set  No.  1 whose  geometric  characteristics  are  described 
in  Section  II-A.  This  rotor  set  had  already  been  tested  as  an  isolated  rotor 
(Reference  1)  so  that  its  performance  in  isolation  was  available  for  purposes 
of  comparison. 

The  results  of  the  investigation  of  the  effect  of  circumferential 
inlet  distortion  on  rotating  stall  are  presented  in  the  following  paragraphs. 

The  calibration  of  the  stationary  multi-lobe  distortion  screens  is  presented 
first.  This  is  followed  by  a presentation  of  the  results  obtained  with  sta- 
tionary distortion.  Finally  the  results  of  the  rotating  distortion  studies 
are  presented. 

1.  Calibration  of  Stationary  Distortion  Screens 
— — —.1—  . — — — . . - 

The  notation  used  m the  distortion  studies  is  shown  in  Figure  in 
and  a sketch  of  the  annular  cascade  with  the  appropriate  geometric  dimensions 
is  given  in  Figure  11.  As  discussed  in  Section  1 1 - A - 2 and  illustrated  ir 
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Figure  11,  the  distortion  screens  were  mounted  on  the  upstream  rotor  hub  so 
that  they  could  be  rotated  about  the  compressor  axis.  In  the  tests  with  sta- 
tionary inlet  distortion,  two  circumferential  distortion  screens  were  used, 
one  with  two  lobes  and  one  with  four  lobes.  These  screens  are  sketched  in 
Figure  12.  Figure  12(a)  shows  the  symmetric  two-lobe  pattern  (Configuration  A) 
Each  lobe  has  a circumferential  extent  encompassing  60  degrees.  Figure  12(b) 
shows  the  symmetric  four-lobe  pattern  (Configuration  B) . Each  lobe  of  this 
pattern  encompasses  30  degrees.  Tr  1 pressure  surveys  were  made  in  the 
wakes  of  the  screens  with  a small  total  pressure  rake  which  can  be  traversed 
circumferentially  over  a limited  range  (approximately  30  degrees),  in  combina- 
tion with  the  various  screen  locations.  The  screens  could  be  indexed  to  any 
circumferential  position.  The  data  were  then  analyzed  as  though  the  screens 
were  held  fixed  at  a reference  position  and  the  rake  traversed  over  the  com- 
plete angular  range  of  interest.  Figures  12(a)  and  12(b)  show  the  reference 
locations  of  the  screen  elements. 

The  results  obtained  from  the  total  pressure  surveys  downstream  of 
screen  configuration  A are  shown  in  Figures  13(a)  and  13(b)  for  two  different 
rotor  speeds,  500  and  950  rpm  respectively.  Similar  data  are  shown  in 
Figures  14(a)  and  14(b)  for  screen  configuration  B.  In  all  of  these  figures, 
the  total  pressure  coefficient,  Cp ^ (r,d)  is  defined  as 

H - W 

-,-f 

-r-  d, 

where  H is  the  total  pressure  upstream  of  the  screen 

Ht  is  the  total  pressure  in  the  survey  plane  downstream  of  the 

screen  (but  upstream  of  the  rotor) 

Ur  is  the  mean  axial  velocity  measured  upstream  of  the  screen 

,C  density 


Each  figure  contains  circumferential,  (0),  surveys  for  various  radial  locations 
.ir  . The  origin  for  the  vertical  total  pressure  coordinate  is  displaced  for 
each  in  proportion  to  the  value  of  . 
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The  two  rotor  rpm's  were  investigated  with  each  screen  to  determine 
if  there  was  any  significant  upstream  influence  from  the  rotor  on  the  screen 
generated  distortion  patterns.  At  500  rpm  the  rotor  was  practically  unloaded 
and  at  950  rpm  the  rotor  was  loaded  nearly  to  the  inception  point  for  rotating 
stall.  Comparison  of  Figures  13(a)  and  13(b)  for  configuration  A and  Figures 
14(a)  and  14(b)  for  configuration  B indicates  that  any  upstream  effect  of  the 
rotor  angular  velocity  cn  the  total  pressure  was  insignificant. 

Ir.  general,  the  distortion  screens  behaved  as  expected,  generating 
a relatively  sharp  edged  wake  with  a total  pressure  coefficient  relatively 
close  to  unity.  In  all  cases,  the  hub  and  tip  wall  boundary  layers  are  evident 
in  the  data  for  the  smallest  and  the  largest  values  of  Ar  , where  the  total 
pressure  coefficients  remain  below  zero  even  outside  of  the  screen  wakes. 

Following  the  total  pressure  calibrations,  the  two  distortion  screens 
were  calibrated  for  swirl  angle.  All  swirl  angle  measurements  were  made  with 
a crossed  hot-film  probe  operated  by  two  Thermo- Systems  hot-wire  anemometers. 

The  outputs  from  the  two  anemometers  were  linearized  and  then  integrated  for 
a period  of  sixty  seconds  in  order  to  obtain  an  accurate  time  average  of  the 
flow  angles.  Air  temperature  variations  were  corrected  for  in  the  initial 
calibration  procedure  prior  to  a test  run.  Since  each  data  point  in  the  swirl 
angle  surveys  required  a minimum  elapsed  time  greater  than  sixty  seconds,  the 
swirl  angle  measurements  are  more  limited  in  the  r umber  of  data  points  gathered 
than  were  the  total  pressure  data. 

The  results  of  swirl  angle  surveys  between  the  rotor  and  the  two- lobe 
distortion  screen  (Configuration  A)  are  shown  in  Figures  15(a)  and  15(b)  for 
rotor  speeds  of  500  and  950  rpm  respectively.  Similar  results  are  shown  in 
Figures  16(a)  and  16(b)  for  the  four-lobe  distortion  screen  (Configuration  B) . 
Each  figure  contains  circumferential  (9)  surveys  for  five  different  radial 
locations  (AO-  The  sign  convention  for  swirl  angle  (/3)  is  shown  in  Figure  10. 

The  circumferential  distributions  of  swirl  angle  generated  by  both 
distortion  screens  are  asymmetric,  with  the  data  measured  near  the  hub  showing 
the  largest  variations.  Moreover,  unlike  the  total  pressure  distributions 
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which  were  insensitive  to  rotor  speed,  the  swirl  angle  distortion  increases 
significantly  as  rotor  speed  is  increased.  Almost  all  of  the  increased  dis- 
tortion in  swirl  angle  with  increasing  rotor  speed  occurs  over  the  lower  half 
of  the  screens  (the  lower  edge  of  the  screen  is  the  one  with  the  smallest 
value  of  Q ).  The  lar.  5 variations  in  swirl  angle,  between  the  distortion 
screens  and  rotor  are  similar  to  those  observed  under  similar  conditions  in 
other  experiments.  For  example,  see  Reference  7 which  comments  on  the  results 
obtained  in  several  experiments.  Apparently  the  swirl  angle  variations  are 
predicted  reasonably  well  by  actuator  disc  theory  and  the  variations  decrease 
with  distance  upstream  of  the  rotor. 

2.  Experiments  With  Stationary  Inlet  Distortion 

Following  the  calibration  of  the  two  stationary  distortion  screens, 
the  effect  of  these  inlet  distortions  on  the  steady  state  turning  and  loss 
performance  of  the  rotor  and  their  effect  on  rotating  stall  was  studied.  The 

turning  and  loss  performance  of  the  rotor  is  presented  first. 

Total  pressure  surveys  downstream  of  the  rotor  were  made  with  both 
screen  configurations  for  a wide  range  of  rotor  rpm  and  with  the  mean  axial 
velocity,  U,  , held  constant  at  60  fps.  The  rotor  used  in  these  tests  is 
Rotor  Set  No.  1 of  Reference  1 set  at  a stagger  angle,  6RM  of  40  degrees. 

With  screen  Configuration  A,  the  rotor  rpm's  used  for  the  tests  were  500,  700, 

850,  950,  1000,  1050,  and  1150.  Rotating  stall  occurred  during  the  tests  at 

1050  and  1150  rpm.  With  screen  Configuration  B,  the  rotor  rpm's  were  500, 

700,  850,  950,  1000,  1050,  1150,  and  1200  with  rotating  stall  present  at  the 
three  highest  values.  These  surveys  were  made  with  a total  pressure  rake 
which  spanned  approximately  one-half  of  the  depth  of  the  annulus.  The  com- 
plete annulus  was  spanned  by  taking  two  sets  of  data,  one  with  the  rake  set 
close  to  the  hub,  and  one  with  the  rake  set  close  to  the  tip.  The  rake  in- 
corporated a Conrad  arrowhead  style  yaw  meter  at  its  mid-span.  This  was  used 
to  align  the  rake  with  the  local  flow  direction.  Thus  if  there  were  signif- 
icant deviations  in  flow  angle  between  the  outer  and  inner  portions  of  the 
annulus,  the  two  data  runs  would  allow  approximate  alignment  of  the  rake  to 
account  for  these  deviations. 


The  results  of  the  total  pressure  surveys  are  shown  in  Figures  17(a) 
through  17(h)  for  the  two-lobe  screen  (Configuration  A)  and  in  Figures  18(a) 
through  18(h)  for  the  four-lobe  screen  (Configuration  B) . For  a given  screen 
configuration,  each  figure  corresponds  to  a different  rotor  rpm  and  the  range 
of  rpm's  extends  from  nearly  unloaded  (500  rpm)  to  loading  beyond  the  inception 
point  of  rotating  stall  (above  1000  rpm).  The  total  pressure  coefficient, 

Cpr  ( r,  6)  is  defined  as 


where  Ht  is  the  total  pressure  upstream  of  the  screen 
Hi  is  the  total  pressure  downstream  of  the  rotor 

U,  is  the  mean  axial  velocity  measured  upstream  of  the  screen 


c is  the  air  density 


Each  figure  contains  circumferential  (&)  surveys  for  various  radial 
locations, a r . The  horizontal  reference  line  for  the  total  pressure  coordi- 
nate is  displaced  for  each  value  of  JSr  in  proportion  tc  the  value  of  . 

Note  that  these  reference  lines  do  not  always  correspond  to  Cpr  = 0.  At 
the  higher  rotor  speeds  for  each  configuration,  CP  - 1 was  chosen  as  the 
reference  line  in  order  to  maintain  the  vertical  scale  ratio  for  Jr  in  these 
figures . 


The  complete  total  pressure  surveys  presented  in  Figures  1“  and  18 
show  several  interesting  features.  First,  the  distortion  patterns  generated 
by  the  screens  are  not  completelv  diffused  by  the  rotcr.  The  wakes  are  most 
evident  near  the  hub  (small  ) and  tend  to  be  more  diffused  after  rotating 
stall  inception  (above  1000  rpm).  Second,  the  data  show  distinct  positive 
peaks  in  the  total  pressure  distributions  near  mid-annulus  at  a circumferential 
location,  9 , approximately  ten  degrees  less  than  the  location  of  the  lower 
edge  of  the  screen  (the  lower  edge  is  the  one  with  the  smallest  value  of  6 ) . 
The  rotor  spins  in  the  direction  of  decreasing  9 and  any  deflection  by  the 
rotor  of  the  screen  wakes  should  also  occur  towards  smaller  values  of  9 . 


At  the  smallest  rotor  speeds,  there  is  evidence  of  a corresponding  negative 
peak  in  the  total  pressure  near  the  upper  edge  of  the  screen,  but  this  is 
small  and  dies  out  quickly  with  increasing  rpm.  Moreover,  the  positive  peak 
in  the  total  pressure  near  the  lower  edge  of  the  screen  dies  out  after  ro- 
tating stall  inception  for  both  distortion  screens.  It  is  possible  that  the 
peaks  observed  in  the  total  pressure  distributions  ar*.  a result  of  the  response 
of  the  rotor  to  the  variations  in  inlet  swirl  angle  generated  by  the  distor- 
tion screens.  For  example,  the  large  negative  values  of  inlet  swirl  angle 
shown  in  Figures  IS  and  16  near  the  lower  edges  of  the  screens  nay  be  the 
cause  of  the  positive  peaks  in  downstream  total  pressure  in  this  same  region. 
The  negative  swirl  angles  in  this  region  increase  the  local  angle  of  attack 
of  the  rotor  and  should  increase  the  work  performed  by  the  rotor  and  thus  the 
downstream  total  pressure.  The  converse  is  true  near  the  upper  edges  of  the 
screen  but  the  effect  is  smaller. 

Swirl  angle  distributions  were  measured  downstream  of  the  rotor  only 
for  the  two-lobe  distortion  screen  (Configuration  A).  As  will  be  discussed, 
it  was  found  that  there  were  negligible  differences  in  rotating  stall  incep- 
tion and  properties  after  inception  for  the  two  distortion  screen  configura- 
tions. This  result  justified  limiting  the  completion  of  these  tests  to  one 
screen  configuration.  Circumferential  surveys  were  made  at  five  radial  lo- 
cations for  each  of  four  different  rotor  speeds  (500,  850,  950  and  1050  rpm). 
This  range  of  rotor  speeds  extends  from  nearly  unloaded  (500  rpm)  to  a loading 
beyond  the  inception  point  of  rotating  stall  (above  1000  rpm).  The  measurement 
of  these  swirl  angles  wes  accomplished  with  the  hot-film  instrumentation 
described  previously. 

The  results  of  the  swirl  angle  surveys  are  presented  in  Figures  19(a) 
through  19(d)  with  each  figure  corresponding  to  a different  rotor  speed.  Note 
that  the  horizontal  reference  lines  for  the  swirl  angle  coordinate  do  not 
correspond  to  zero  values  of  swirl  angle,  /3a  , and  that  the  reference  values 
vary  from  figure  to  figure.  The  variable  reference  values  were  selected  in 
order  to  maintain  the  same  vertical  scale  ratio  for  f&  in  these  figures. 
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The  swirl  angles  downstream  of  the  rotor  vary  with  circumferential 
location.  The  largest  deviations  from  the  mean  value  occur  just  before  ro- 
tating stall  inception  (9S0  rpm)  and  decrease  after  inception  (1050  rpm) . 

The  extremes  of  the  circumferential  variations  in  downstream  swirl  angle  are 
approximately  half  of  .-.hose  observed  upstream  of  the  rotor  and  tney  do  not 
appear  to  form  a consistent  trend  which  can  be  associated  with  the  presence 
of  the  distortion  screen  upstream  of  the  rotor. 


For  each  radial  location  at  a given  rotor  speed,  circumferential  j 

: averages  of  have  been  calculated.  Since  the  data  only  extend  between 

24  deg  s?  9 6 168  deg,  the  full  0 to  180  degree  circumferential  range  required  < 

i 1 

for  computing  circumferential  averages  (the  distortion  screen  is  a symmetric  j 

two-lobe  pattern)  was  obtained  by  assuming  a linear  variation  in  /3z  for  the 

region  where  5X  was  not  measured.  The  resulting  circumferential  average  j 

values  of  the  swirl  angle,  /3X  , are  compared  in  Figure  20  with  similar  data 

for  the  isolated  rotor  (Ref.  1).  In  Figure  20,  two  degrees  have  been  added  ] 

to  the  data  measured  with  the  screen  in  place.  This  has  been  done  in  order 

to  allow  direct  comparison  with  the  shapes  of  the  radial  distributions  for  the 

isolated  rotor.  Isolated  rotor  data  are  not  available  for  a rotor  speed  of 

850  rpm.  Thus  a direct  comparison  for  this  speed  cannot  be  made.  However, 

the  isolated  rotor  data  for  rotor  speeds  of  800  and  900  rpm  are  shown  in  the 

f figure  to  provide  an  indirect  comparison.  At  1050  rpm,  both  the  isolated 

| rotor  and  the  rotor  with  inlet  distortion  were  undergoing  rotating  stall. 

! ; 

It  is  evident  from  inspection  of  Figure  20  that  the  radial  distribu- 
tions of  fi  with  and  without  the  distortion  screen  are  the  same  shapes  to 
within  the  probable  experimental  accuracy  of  the  data.  In  addition,  it  would 
appear  that  the  radial  distributions  of  /3  with  the  screen  in  place  are  an 
almost  constant  two  degrees  less  than  those  observed  on  the  isolated  rotor. 

| In  general,  the  effect  of  the  distortion  screen  on  the  flow  turning  performance 

[ of  the  rotor  is  very  small.  The  rotating  stall  stability  theory  in  Section  III 

requires  an  overall  average,  /3  , obtained  by  integrating  the  circumferential 

averages  /$x  in  the  radial  direction.  In  view  of  the  consrant  two  degree  dif- 
ferences in  the  circumfereintial  averages,  it  is  sufficiently  accurate  in  this 
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case  to  assume  that  the  isolated  rotor  measurements  diminished  by  a constant 
two  degrees  can  be  used  for  input  to  the  theory. 

It  is  somewhat  moie  difficult  to  compare  the  loss  performance  of  the 
rotor  with  and  without  inlet  distortion.  In  principle,  such  a comparison 
requires  that  the  radial  distributions  of  total  pressure  measured  in  the 
absolute  coordinate  system  (as  in  Figures  17  and  18)  be  converted  to  a co- 
ordinate system  relative  to  the  rotor  and  then  integrated  radially  as  in 
Reference  1.  The  swirl  velocities  downstream  of  the  rotor  enter  into  this 
integration  procedure.  However  in  view  of  the  small  constant  difference  in 
the  radial  distribution  of  swirl  angle,  p1  , observed  between  the  isolated 
rotor  and  the  rotor  with  inlet  distortion,  a valid  comparison  of  rotor  loss 
performance  can  be  made  on  the  basis  of  the  total  pressure  measurements  made 
in  the  absolute  coordinate  system.  This  has  been  done  and  the  results  are 
presented  in  Figure  21. 

In  Figure  21,  the  overall  total  pressure  difference  across  the  rotor 
is  plotted  in  coefficient  form  as  a function  of  flow  coefficient.  Both  the 
data  for  the  isolated  rotor  and  for  the  flows  with  inlet  distortion  are  shown. 
The  notation  in  the  figure  is  as  follows: 


where: 


q.t  ■ Dynamic  pressure  corresponding  to  average  axial  velocity, ( ~ U*) 

• Dynamic  pressure  corresponding  to  average  velocity  relative 
to  the  rotor,  (Ut*  + ) 
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» Area  Weighted  Average  of  Total  Pressure  at  Rotor  Face 
(Measured  in  absolute  coordinate  system) 

M 

PT  • Area  Weighted  Average  of  Total  Pressure  Leaving  Rotor 
(Measured  in  absolute  coordinate  system) 

U0  ■ Area  Weighted  Average  of  Axial  Velocity  Upstream  of  Rotor 

Wb  ■ Area  Weighted  Average  of  Rotor  Velocity. 

The  average  total  pressures  were  computed  by  numerically  integrating 
the  measured  total  pressure  data  presented  in  Figures  13,  14,  17  and  18  for 
the  cases  with  inlet  distortion.  These  data  include  the  effect  of  the  boundary 
layers  developed  on  the  hub  and  outer  casing  of  the  annular  cascade  both  up- 
stream and  downstream  of  the  rotor.  The  total  pressure  data  for  the  isolated 
rotor  (Reference  1)  included  the  wall  boundary  layers  downstream  of  the  rotor 
but  not  upstream  of  the  rotor.  Tn  Figure  21,  these  original  data  for  the 
isolated  rotor  were  corrected  to  include  the  upstream  boundary  layers  by  as- 
suming that  the  dimensionless  total  pressure  data  measured  between  the  dis- 
tortion screen  lobes  were  representative  of  the  upstream  boundary  layers  on 
the  isolated  rotor.  This  procedure  should  be  fairiy  accurate  since  the  geometry 
of  the  annular  cascade  was  identical  except  for  the  segments  of  the  distortion 
screens.  Moreover  the  data  from  both  screen  configurations  provided  the  same 
value  for  the  boundary  layer  correction  to  the  isolated  rotor. 

It  can  be  seen  from  figure  21  that  the  dimensionless  total  pressure 
differences  across  the  rotor  are  nearly  identical  for  all  of  the  data  except 
in  the  regions  where  rotating  stall  is  occurring.  For  application  of  the  ro- 
tating stall  theory,  it  is  only  permissible  to  use  the  data  measured  prior 
to  rotating  stall  inception.  These  regions  are  labelled  "no  rotating  stall" 
in  Figure  2i  and  the  small  differences  in  the  data  are  within  the  range  of 
the  probable  experimental  scatter.  Hence  it  is  concluded  that  prior  to  ro- 
tating stall  inception  the  average  loss  performance  of  the  rotor  with  inlet 
distortion  was  essentially  the  same  as  that  found  on  the  isolated  rotor. 
Moreover,  as  can  be  seen  in  Figure  21,  rotating  stall  inception  occurred  at 


almost  identical  values  of  flow  coefficient,  , on  the  isolated  rotor  and 
on  the  rotor  with  circumferential  inlet  distortion. 

In  addition  to  rotating  stall  inception  points,  the  numbers  of 
cells  and  propagation  velocities  of  rotating  stall  were  measured  for  both  dis- 
tortion screens.  These  data  were  obtained  through  the  use  i f two  static 
pressure  taps  on  the  outer  casing  at  the  approximate  axial  location  of  the 
rotor  quarter- chord.  The  two  static  taps  were  separated  by  seven  rotor  blade 
spaces  in  the  circumferential  direction.  The  rotor  contains  46  blades;  thus 
the  circumferential  separation,  A 9 , was  54.8  degrees.  With  the  two-lobe 
distortion  screei,  data  were  taken  for  two  different  screen  circumferential 
locations  with  rospect  to  the  static  pressure  taps.  This  was  done  in  order 
to  determine  if  the  location  of  the  pressure  sensors  relative  to  the  distor- 
tion screen  has  an  effect  on  the  detection  of  rotating  stall  inception  and 
its  properties  after  inception.  The  location  of  the  pressure  sensors  with 
respect  to  the  distortion  screen  had  no  detectable  effect  on  the  inception  or 
properties  of  rotating  stall.  However,  at  rotor  rpm's  well  below  rotating 
stall,  the  pressure  fluctuations  due  to  blade  passage  were  affected  by  the 
relative  location  of  the  pressure  sensors.  A pressure  sensor  in  the  wake  of 
the  distortion  screen  showed  larger  blade  passage  pressure  f’uctuations  than 
a sensor  outside  of  the  screen  wake.  As  rotor  speed  was  increased  toward 
rotating  stall  inception,  the  blade  passage  fluctuations  outside  of  the  screen 
wake  increased  more  rapidly  than  those  inside  the  wake  until,  at  inception, 
the  blade  passage  pressure  fluctuations  were  approximately  equal  inside  ana 
outside  of  the  screen  wake.  Moreover,  rotating  stall  pressure  fluctuations 
were  also  approximately  equal  inside  and  outside  of  the  screen  wake. 

The  measured  rotating  stall  propagation  velocities  are  shown  in 
Figure  22  for  both  distortion  screens.  Previously  obtained  data  (Reference  1) 
for  the  rotor  with  no  inlet  distortion  are  shown  for  comparison.  Note  that 
the  propagation  velocities,  Vp  are  presented  in  a coordinate  system  fixed  to 
the  rotor  blades  and  also  that  the  propagation  velocities  have  been  made 
dimensionless  by  the  mean  inlet  swirl  velocity,  , relative  to  the  rotor. 


28 


The  absolute  rotating  stall  propagation  velocity  is  related  to  the 
relative  velocity  by  = l-(VP/rtp)  where  \/p  is  the  absolute  propaga- 

tion velocity  and  Wp  is  the  rotor  velocity.  The  horizontal  coordinate  of 

a 

Figure  22  is  the  overall  inlet  swirl  angle,  /S,  , relative  to  the  rotor.  The 
value  of  3,m  was  calculated  as  in  Reference  1 (Equation  10)  and  assumes  that 
local  deviations  in  the  relative  inlet  swirl  angle  due  to  the  presence  of  the 
distortion  screens  integrate  to  zero  around  the  annulus.  The  number  of  stall 
cells  are  indicated  in  Figure  22  by  the  numbers  near  the  data  points. 

The  data  presented  in  Figure  22  show  only  minor  differences  in  the 
inception  and  properties  of  rotating  stall  with  and  without  inlet  distortion. 
In  all  cases,  rotating  stall  inception  was  slightly  hysteretic,  that  is  as 
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rotor  rpm  was  increased  ( 3t  increasing)  rotating  stall  occurred  at  a hi,  r 
value  of  rpm  than  the  value  at  which  it  disappeared  on  decreasing  the  rotor 

rpm  ( 3 decreasing).  With  both  distortion  screen  configurations,  inception 

* ^ 

occurred  at  the  same  value  of  ,3,  and  this  value  was  less  than  one  degree 

_ <* 

larger  than  the  value  of  ,3  at  inception  with  no  distortion  screen.  The 
numbers  of  stall  cells  observed  were  similar  in  all  cases,  with  the  changeover 
from  one  to  two  cells  occurring  at  identical  values  of  3.  . With  the  dis- 

tortion screens  in  place,  relative  propagation  velocities  were  slightly  higher 
than  those  observed  with  no  distortion  but  the  overall  scatter  in  the  data  is 
larger  than  the  apparent  di  f fer '->ces  i.i  propagation  velocity.  In  general  it 
can  be  stated  that  the  presence  of  the  distortion  screens  had  a negligible 
effect  on  the  inception  of  rotating  stall  and  on  its  properties  after  incep- 
tion. 


In  summary,  the  comparisons  of  steady  state  turning  and  loss  per- 
formance and  of  rotating  stall  inception  and  properties  on  a rotor  with  and 
without  circumferential  inlet  distortion  has  shown  only  very  small  differences 
in  all  quantities  The  largest  difference  is  a small  constant  difference  in 
flow  turning  performance  between  the  isolated  rotor  and  the  rotor  with  inlet 
distortion.  The  nearly  identical  steady  state  turning  and  loss  performance 
of  the  rotor  with  and  without  inlet  distortion  may  be  fortuitous.  However, 
given  these  steady  state  performance  similarities,  it  is  worth  noting  that 
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rotating  stall  inception  and  its  propei*ies  after  inception  were  also  similar 
with  and  without  inlet  distortion.  The  implications  of  there  results  are 
discussed  in  Section  III  where  the  experimental  results  are  compared  to  the- 
oretical predictions. 

3.  Experiments  With  Rotating  Distortion 

The  final  tests  performed  in  these  studies  of  rotating  stall  were 
designed  to  investigate  whether  there  is  a difference  in  rotating  stall  in- 
ception conditions  when  a circumferential  inlet  distortion  pattern  rotates 
at  different  speeds  about  the  compressor  axis.  The  multi-lobe  distortion  in- 
vestigations mentioned  previously  (Reference  3 and  4)  have  suggested  that 
inception  may  be  sensitive  to  a dimensionless  reduced  frequency  based  upon 
the  time  of  passage  of  a rotor  blade  through  the  distortion  pattern.  If  this 
is  true,  it  should  be  possible  to  observe  changes  in  the  inception  point  on 
a given  rotor  subjected  to  a given  inlet  distortion  pattern  which  has  different 
rotational  speeds  about  the  compressor  axis.  The  changes  will  arise  because 
the  passage  time  of  a rotor  blade  through  the  lobes  of  the  distortion  patten 
will  change  as  the  distortion  pattern  changes  speed.  Evidence  obtained  in 
this  way  would  be  subject  to  less  ambiguity  than  that  inferred  from  the  multi- 
lobe distortion  investigations  because  the  circumterenti ai  extent  of  the  dis- 
tortion pattern  will  remain  essentially  unchanged. 

These  tests  were  performed  with  Rotor  Set  No.  1 and  the  two-lobe 
distortion  screen  (Configuration  A)  sketched  in  Figure  12(a).  The  configura- 
tion of  the  annular  cascade  was  the  same  as  that  used  for  the  tests  with 
stationary  distortion  (Figure  11).  The  procedure  used  in  these  tests  was  to 
hold  the  rotor  speed  constant  at  1000  rpm  and  study  rotating  stall  inception 
at  a series  of  distortion  screen  rotational  speeds.  At  each  distortion  screen 
speed,  inception  was  incurred  by  gradually  reducing  the  mass  flow  through  the 
rotor  from  a high  value  to  a low  value. 

Several  methods  for  detecting  the  presence  of  rotating  stall  were 
investigated.  The  most  informative  data  were  obtained  by  recording  the  unsteady 
total  pressure  fluctuations  at  mid-annulus  immediately  upstream  of  the  rotor. 
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These  records  display  the  distortion  screen  wakes  with  the  fluctu- 
ations due  to  rotating  stall  superimposed  upon  them.  It  is  these  total 
pressure  records  which  are  used  in  the  following  presentation. 

The  experimental  results  are  presented  as  a series  of  Brush  Recorder 
strips  in  Figures  23  through  27.  Each  individual  record  is  composed  of  two 
recorder  strips.  The  upper  strip  shows  the  dynamic  pressure,  , measured 
in  the  constant  area  annulus  well  upstream  of  the  distortion  screen.  (Sta- 
tion 0 in  Figure  11).  The  mean  axial  velocity,  U0  , through  the  test  section 
is  proportional  to  the  square  root  of  . The  value  of  tf  , at  rotating  stall 
inception  in  combination  with  the  constant  rotor  speed  (1000  rpm)  is  suffi- 
cient to  define  the  flow  coefficient,  t?  , at  inception.  The  lower  strip 
shows  the  total  pressure  difference,  A PT  between  two  probes,  one  mounted  at 
mid-annulus  well  upstream  of  the  distortion  screen  (Station  0 in  Figure  11) 
and  one  mounted  at  mid-annulus  just  upstream  of  the  rotor  (Station  1 in 
Figure  11).  The  connection  between  the  pressure  transducer  and  the  total 
pressure  probe  at  Station  1 was  kept  as  short  as  possible  to  obtain  a good 
frequency  response.  The  frequency  response  was  not  measured  but  the  results 
indicate  that  it  was  sufficient  to  display  the  distortion  screen  wakes  and 
fluctuations  due  to  rotating  stall. 

In  each  of  Figures  23  through  27,  the  vertical  scale  factors  for 
and  APt  are  shown  on  the  right  side  of  the  appropriate  record  strips. 

The  scale  factor  for  c\.t  is  the  same  in  all  records  (0.4  to  1.4  inches  of 
water).  Two  different  scale  factors  were  used  for  the  total  pressure  dif- 
ference, A PT  . The  first  experiments,  which  were  performed  with  the  rotor 
and  screen  rotating  in  opposite  directions,  used  a full  scale  range  from  0 to 
-2  inches  of  water.  (Negative  values  indicate  a total  pressure  drop  across 
the  distortion  screen.)  In  later  tests  with  the  rotor  and  screen  rotating  in 
the  same  direction,  this  range  was  found  to  be  too  small.  Thus  the  full 
scale  range  for  A PT  was  changed  to  0 to  -5  inches  of  water.  This  change  in 
scale  should  be  kept  in  mind  when  inspecting  the  A PT  records.  In  all  cases, 
time  is  increasing  from  left  to  right  and  its  magnitude  is  given  between  the 
and  &PT  record  strips.  Changes  in  time  s^ale  during  a record  are  also 
indicated. 
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The  first  sec  of  records,  shown  in  Figure  23,  illustrate  rotating 
stall  inception  with  the  distortion  screen  stationary.  These  records  cor- 
respond to  the  test  conditions  used  in  the  stationary  distortion  screen  ex- 
periments described  in  Section  II-C-2.  The  measurements  made  during  the  sta 
tionary  screen  experiments  were  more  comprehensive  and  used  rotating  stall 
detection  techniques  which  have  been  proven  in  past  studies.  Thus  it  is  of 
interest  to  compare  results  obtained  from  Figure  23  with  the  results  of  the 
stationary  distortion  screen  experiments.  This  is  done  in  the  following  dis 
cuss ion. 


Figure  23(a)  was  obtained  with  the  downstream  total  pressure  probe 
(Station  1)  in  the  free-stream  area  between  the  two  distortion  screen  lobes. 
Figure  23(b)  was  obtained  with  the  same  ^otal  pressure  probe  in  the  wake  of 
one  of  the  distortion  screen  lobes.  The  presence  of  rotating  stall  is  obvi- 
ous in  both  the  <^t  and  the  APT  record  strips.  The  steady  total  pressure 
drop  in  the  distortion  screen  wake  prior  to  rotating  stall  inception  can  be 
seen  on  the  far  left  side  of  Figure  23(b).  The  flow  coefficient  <J>  define^ 
in  Section  II-A-2  can  be  calculated  from  the  rotor  speed  and  the  value  of 
at  inception.  These  calculated  values  of  df  are  0.509  and  0.496  for  Figures 
23(a)  and  23(b)  respectively.  These  values  compare  well  with  che  stationary 
distortion  data  shown  in  Figure  21. 

The  expanded  time  portions  of  the^PT  records  on  the  right  side  of 
Figure  23  show  that  with  stationary  distortion  rotating  stall  is  highly  regu- 
lar with  a frequency  of  8 Hz.  This  corresponds  to  one  stall  cell  moving  at 
an  absolute  propagation  velocity  of  48  percent  of  the  rotor  speed  or  a rela- 
tive velocity  of  52  percent  of  the  rotor  speed.  The  relative  value  falls 
within  the  range  of  values  observed  in  the  earlier  tests  with  stationary  dis- 
tortion. (Figure  22,  Configuration  A).  The  main  difference  between  the  total 
pressure  records  in  Figures  23(a)  and  23(b)  is  that  the  distortion  screen 
wake  is  evident  both  before  and  after  rotating  stall  inception  in  Figure  23(b). 
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In  view  of  the  good  agreeme.  between  the  results  for  rotating  stall 
inception  inferred  from  the  data  presented  in  Figure  23  and  the  results  ob- 
tained in  the  more  comprehensive  measurements  presented  in  Section  II-A-2  for 
stationary  distortion,  it  is  concluded  that  the  technique  used  to  generate 
Figure  23  can  be  used  to  study  the  effect  of  rotating  distortion  on  inception 
of  rotating  stall.  The  results  obtained  with  rotating  distortion  are  pre- 
sented in  the  following  paragraphs. 

Figure  24  presents  a series  of  Brush  record  strips  taken  with  the 
distortion  screen  rotating  in  the  direction  opposite  that  of  the  rotor. 

Parts  a,  b and  c of  Figure  24  show  the  results  obtained  at  three  different 
rotational  speeds  of  the  distortion  screen.  The  full  series  of  tests  covered 
the  range  of  screen  rpm's  between  0 and  -1000  in  increments  of  100  rpra.  As 
with  the  stationary  screen,  inception  of  rotating  stall  is  obvious  in  both 
the  APr  and  the  q.t  record  strips. 

The  left  side  of  Figure  23  shows  conditions  prior  to  inception. 

The  A PT  records  indicate  that  the  magnitude  of  the  defect  in  the  wakes  of  the 
rotating  screen  increased  with  increasing  screen  rpm.  This  effect  will  be 
illustrated  more  clearly  in  later  figures  where  the  average  value  of  was 
held  constant.  In  this  figure  was  varied  in  order  to  induce  inception. 

The  cf0  records  prior  to  rotating  stall  also  show  fluctuations  at  the  passage 
frequency  of  the  distortion  screen  lobes.  This  indicates  that  the  screen 
generated  disturbances  extend  well  upstream  of  the  rotating  screen.  The  mag- 
nitude of  these  fluctuations  are  probably  larger  than  those  in  the  record 
strips  since  the  qa  transducer  was  connected  to  the  upstream  total  and  static 
pressure  probes  by  long  tubing.  For  the  purpose  of  these  experiments,  only 
the  time  average  value  of  <].g  is  required  at  rotating  stall  inception.  This 
was  obtained  by  numerically  averaging  the  fluctuating  ^ records  at  inception. 

Inspection  of  the  expanded  time  scale  records  on  the  right  of 
Figure  24  suggests  that  rotating  stall  is  almost  completely  independent  of 
the  distortion  screen  disturbances.  This  is  particularly  apparent  at  the 
two  lower  screen  ri  's  where  the  ^PT  signals  from  rotating  stall  are  easily 


distinguishable  from  the  screer  wakes  because  the  stall  signals  are  substan- 
tially larger  than  the  wake  signals.  The  rotating  stall  signals  show  the 
same  characteristics  as  those  obtained  with  the  stationary  screen  (Figure  23). 
The  screen  wake  signals  are  simple  superimposed  on  the  rotating  stall  signals. 
Moreover,  the  fluctuating  ^ records  show  that  the  stall  frequency  is  nearly 
independent  of  screen  rpm.  Thi,  frequency  decreased  a very  small  amount  with 
increased  screen  rpm,  fro-.  S He  at  0 screen  rpm  (0  Hi)  to  5.8  Hz  at  -800  screen 
rpm  (26.7  Hz). 

Inspection  of  the  inlet  dynamic  pressure,  cf  , records  at  rotating 
stall  inception  shows  that  as  negative  screen  rotation  increased,  the  average 
value  of  if'  , and  hence  the  average  flow  coefficient  if  , increased  by  a small 
amount.  The  quantitative  effect  of  screen  rotation  on  the  average  flow  co- 
efficient at  inception  will  be  presented  shortly  along  with  results  obtained 
with  positive  screen  rotation.  At  this  point  it  is  sufficient  to  note  that 
the  rotating  stall  inception  points  defined  by  and  also  the  frequencies 
of  rotating  stall  after  inception  were  almost  independent  of  distortion 
screen  rotation  in  a direction  opposi'e  that  of  the  rotor. 

Sample  records  obtained  with  the  distortion  screen  rotating  in  the 
same  direction  as  the  rotor  are  shown  in  Figure  25.  Note  that  the  L P„  records 
have  a different  full  scale  range  than  those  shown  in  Figure  2d.  As  with  the 
negative  screen  rotation,  the  full  series  of  tests  covered  the  complete  range 
of  screen  rpre's  between  - and  1^00  m increments  of  100  rpm.  Additional  re- 
sults were  recorded  at  50  and  550  rpra  on  the  screen  because  of  some  interest- 
ing features  in  the  overa.i  results. 

Prior  to  inception  of  rotating  stall,  the  ^ and  £ Pr  records  in 
'igure  25  show  screen  generatud  fluctuations  similar  to  those  observed  at 
negative  screen  rotations.  However  the  average  value  of  at  inception 
varies  considerably  with  positive  screen  rotation.  The  complete  set  of 
records  taken  with  positive  screen  rotation  showed  that  there  is  a large  peak 
in  the  inception  values  of  ^ . This  peak  is  centered  about  a screen  speed 
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of  approximately  5S0  rpm.  At  lower  and  higher  screen  speeds  the  inception 
values  of  ^ decrease.  This  result  is  different  from  that  obtained  with 
negative  screen  rotational  speeds  where  the  value  of  at  inception  was 
almost  independent  of  screen  speed. 

After  stall  inception,  the  APr  records  display  a variable  behavior. 
At  low  screen  speeds  (Figure  25(a))  the  rotating  screen  wakes  appear  to  inter- 
fere destructively  with  rotating  stall,  limiting  the  occurrence  of  large 
amplitude  stall  to  intermittent  intervals  in  time.  This  destructive  inter- 
ference appeared  to  reach  a maximum  at  a screen  speed  of  approximately 
300  rpm  (Not  shown).  At  screen  rotational  speeds  near  500  rpra,  (Figure  25(b)) 
the  screen  generated  disturbances  become  nearly  in  phase  with  rotating  stall 
and  appear  to  reinforce  the  occurrence  of  rotating  stall.  The  stall  frequency 
under  these  conditions  appears  to  double,  corresponding  to  a two  cell  ro- 
tating stall  pattern  appropriate  to  the  two  lobe  rotating  distortion  screen. 

At  screen  speeds  well  above  500  rpm,  destructive  interference  begins  to  occur 
again  (Figure  2S(c)). 

The  screen  speeds  at  which  rotating  stall  is  reinforced  correspond 
to  the  peak  values  of  ^ measured  at  inception  while  the  destructive  inter- 
ference corresponds  to  decreasing  values  of  at  inception.  In  the  regions 
of  destructive  interference,  it  is  possible  to  distinguish  the  occurrence  of 
a small  amplitude  and  a large  amplitude  rotating  stall  at  different  values  of 
. with  the  small  amplitude  inception  occurring  at  higher  values  of  Cf.t  . 
These  inception  points  have  been  labeled  in  Figure  25  in  cases  where  they  are 
distinguishable.  It  is  worth  noting  that  the  results  obtained  with  negative 
screen  rotation  (Figure  24 j showed  only  largo  amplitude  stall  inception, 
while  the  stationary  screen  (Figure  23)  showed  small  amplitude  inception  which 
rapidly  changed  to  large  amplitude  without  a change  in  . 

In  summary,  the  results  presented  in  Figure  24  and  25  for  negative 
and  positive  values  of  distortion  screen  rotational  speeds  suggest  that  ro- 
tating stall  displays  an  effect  similar  to  resonance  with  the  maximum  effect 
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centered  at  a screen  speed  of  approximately  550  rpm  (55  percent  of  the  rotor 
speed) . In  a rotor  fixed  coordinate  system  this  corresponds  to  a speed  ratio 
of  0.45  which  agrees  very  well  with  the  natural  rotating  stall  propagation 
speeds  observed  on  the  isolated  rotor  and  on  the  rotor  with  stationary  dis- 
tortion Configuration  A (See  Figure  22) . 

Some  additional  records  which  display  the  similarity  to  resonance 
are  shown  in  Figures  26  and  27.  In  these  figures  the  average  inlet  dynamic 
pressure  was  held  constant  at  one  inch  of  water  and  a series  of  records  were 
obtained  for  fine  increments  in  screen  rotational  speed.  Rotor  speed  was 
held  constant  at  1000  rpm.  The  result  is  somewhat  similar  to  inspecting  the 
amplitude  versus  frequency  response  of  a forced  resonant  system  while  holding 
the  forcing  amplitude  nearly  constant.  The  correspondence  is  not  exact  since 
the  forcing  function  amplitude,  APt  , increased  with  increasing  screen  speed 
and  the  response  depends  on  the  constant  value  of  selected.  In  this  case, 
was  selected  to  be  above  the  value  required  for  rotating  stall  inception 
except  in  the  region  of  the  resonant  peak. 

The  records  obtained  with  negative  screen  speeds  are  shown  in 
Figure  26.  Rotating  stall  was  not  occurring  in  any  of  these  records.  The 
main  feature  to  be  observed  is  the  increase  in  amplitude  of  the  total  pres- 
sure defect,  APr  , as  screen  speed  is  increased.  This  occurred  even  though 
the  value  of  ^ was  held  constant.  It  is  believed  to  be  a result  of  increas- 
ing swirl  angle  relative  to  the  screen  as  the  screen  speed  is  increased.  The 
screen  is  composed  of  flat  plates  with  holes  drilled  through  them  to  attain 
the  required  porosity.  In  this  type  of  structure,  changing  the  flow  angle 
from  a direction  normal  to  its  face  results  in  increased  losses  due  to 
changes  in  leading  edge  separation  on  the  holes.  The  effect  should  be  inde- 
pendent of  the  direction  of  rotation.  Similar  results  have  been  observed  on 
honey  comb  structures.  In  any  event,  as  noted  previously  the  increase  in 

amplitude  with  negative  screen  speed  did  not  have  much  influence  on  ro- 
tating stall  inception  or  on  its  properties  after  inception. 
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The  records  obtained  with  positive  screen  rotational  speeds  are 
shown  in  Figure  27.  At  screen  rpa's  of  300  and  less  (Figures  27(a),  (b)  and 
(c))  rotating  stall  did  not  occur  while  at  400  rpm  there  is  evidence  of  an 
occasional  snail  amplitude  rotating  stall  in  the  &PT  record  (Figure  27(d)). 

At  500  and  550  rpn  (Figures  27(e)  and  (f)),  large  amplitude  rotating  stall 
occurs  regularly  and  is  nearly  in  phase  with  the  screen  generated  disturbances. 
At  rotational  speeds  above  550  rpn,  rotating  stall  beeches  intermittent  and 
of  smaller  amplitude  until  at  800  rpn,  detectable  rotating  stall  has  disap- 
peared. In  the  absence  of  rotating  stall,  the  screen  generated  disturbances 
in  APr  increase  with  screen  speed  in  approximately  the  same  fashion  as  ob- 
served with  negative  screen  speeds  (Note  the  change  in  full  scale  range  of 
the &Pr  records  between  Figures  26  and  27.)  At  high  positive  screen  speeds, 
the  expanded  time  records  for  aPt  display  a change  in  shape  but  no  significant 
changes  in  amplitude.  These  changes  may  be  due  to  an  interaction  between 
the  distortion  and  the  rotor  as  the  distortion  speed  approaches  the  rotor 
speed.  However,  they  are  defini.oly  not  evidence  of  large  amplitude  rotating 
stall  and  are  not  believed  to  be  evidence  of  small  amplitude  rotating  stall. 

As  noted  previously,  quantitative  values  for  the  average  flow  coef- 
ficient, $ , at  rotating  stall  inception  can  be  calculated  from  the  average 
value  of  q,  at  inception.  This  has  been  done  for  all  of  the  records.  The 
results  are  shown  in  Figure  28  plotted  as  a function  of  distortion  screen  rpm 
divided  by  the  rotor  rpm.  When  presented  in  this  way,  the  resulting  curve 
is  similar  to  that  of  a damped  harmonic  oscillator  with  its  fundamental  fre- 
quency occurring  at  a screen/rotor  velocity  ratio  of  about  0.55.  There  are 
no  data  available  to  define  the  peak  of  the  resonant  curve  because  the  max- 
imum flow  coefficient  attainable  in  the  annular  cascade  was  reached  before 
rotating  stall  disappeared  in  this  region.  The  maximum  attainable  value  of 
$ is  shown  as  a leavy  dashed  line  in  Figure  28.  Points  A and  B in  the  peak 
region  were  obtained  at  the  maximum  attainable  S rather  than  at  inception. 

The  corresponding  notes  for  A and  B are  included  to  describe  the  rotating 
stall  characteristics  at  this  maximum  value  of  . The  data  presented  in 
Figures  26  and  27  for  constant  average  correspond  to  a constant  average 
flow  coefficient,  <$  , of  0.586  in  Figure  28.  Thus  Figures  27(e),  (f)  and  (g) 
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display  large  amplitude  rotating  stall  at  screen/rctor  velocity  ratios  of 
O.S,  0.55,  and  0.6  respectively,  while  Figures  27(d)  and  (h)  displav  small 
amplitude  rotating  stall  for  screen/rotor  velocity  ratios  of  0.4  and  0.7 
respectively. 

In  previous  experiments,  tne  propagation  velocities  of  rotating 
stall  on  the  isolated  rotor  and  on  the  rotor  with  stationary  inlet  distortion 
were  measured.  The  results  obtained  for  these  naturally  occurring  rotating 
stalls  have  been  presented  in  Figure  22  for  a frame  of  reference  fixed  to  the 
rotor.  The  range  of  these  natural  stall  propagation  speeds  are  shown  in 
Figure  28  in  the  stationary  coordinate  system  appropriate  to  this  figure. 

Note  that  the  peak  in  the  rotating  stall  inception  curve  occurs  at  screen/ 
rotor  speed  ratios  which  agree  almost  exactly  with  the  natural  rotating  stall 
propagation  speed. 

The  general  characteristics  for  the  rotating  stall  inception  data 
shown  in  Figure  28  agree  with  the  observations  made  previously  regarding  the 
cnaracteristics  of  rotating  stall  when  it  does  occur.  These  may  be  summarized 
as  follows.  At  negative  and  zero  screen/rotor  velocity  ratios,  rotating  stall 
amplitude  and  frequency  are  almost  completely  independent  of  screen  velocity. 

At  low  positive  screen/rotor  velocity  ratios  there  is  destructive  interference 
between  the  distortion  and  the  rotating  stall.  Rotating  stall  becomes  inter- 
mittent and  occurs  as  a small  amplitude  stall  followed  by  a larger  amplitude 
stall  at  a lower  value  of  d?  . The  maximum  destructive  interference  was  ob- 
se~\ed  with  a screen/rotor  velocity  ratio  of  0.3.  The  samj  value  shows  the 
maximum  delay  in  inception  of  large  amplitude  rotating  stall  (lowest  vali  n 
of  $ ) . Near  the  peak  of  the  inception  curve,  rotating  stall  becomes  large 
in  amplitude  and  approximately  in  phase  with  the  distortion  screen  wakes. 
Moreover  the  rotating  stall  frequency  doubles  from  its  previous  nearly  constant 
value.  This  indicates  that  the  number  of  stall  cells  has  changed  from  one 
to  two  in  correspondence  with  the  number  of  lubes  on  the  distortion  screen. 

At  screen/rotor  velocity  ratios  larger  than  the  value  at  the  peaks  of  the  in- 
ception curve,  rotating  stall  again  becomes  intermittent  and  the  amplitude 
decreases.  Here  again,  large  amplitude  rotating  stall  is  preceded  by  small 
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amplitude  rotating  stall  and  the  flow  coefficient  at  inception  decreases 
with  increasing  distortion  screen  speed. 

The  above  observations  on  destructive  interference  and  reinforce- 
ment and  the  resonant- like  peak  response  of  the  inception  curve  contradict 
the  hypothesis  that  there  is  a minimum  time  required  to  stall  the  rotor 
blades,  and  thus  that  the  time  required  for  a rotor  blade  to  pass  through  a 
distortion  cell  will  have  a significant  influence  on  rotating  stall.  If 
the  hypothesis  were  valid,  rotating  stall  inception  would  occur  at  increasingly 
larger  values  of  flow  coefficient  as  the  screen  speed  approached  the  rotor 
speed.  Instead  the  data  show  a decrease  in  flow  coefficient  near  a screen/ 
rotor  speed  ratio  of  one.  Moreover,  the  destructive  interference  observed 
in  the  data  cannot  be  reconciled  with  the  time  to  stall  hypothesis. 

In  many  ways  the  results  obtained  in  the  rotating  distortion  ex- 
periments are  similar  to  the  results  observed  in  the  damping  experiments  of 
Reference  1.  In  those  experiments  a disturbance  was  generated  by  rotating 
a small  flat  plate  about  a radial  axis.  The  plate  was  situated  between  two 
of  the  blades  in  a stator  row.  It  was  found  that  rotating  stall  is  almost 
completely  independint  of  the  disturbance  excitation  frequency  unless  the 
disturbance  frequency  and  the  frequency  of  naturally  occurring  rotating  stall 
were  closely  matched.  A small  amount  of  destructive  interference  between 
the  forcing  disturbance  and  rotating  stall  was  also  observed  at  disturbance 
generator  frequencies  removed  from  the  natural  rotating  stall  frequency. 

Thus  both  the  previous  damping  experiments  and  the  current  rotating 
distortion  experiments  suggest  that  the  response  of  a blade  row  to  circum- 
ferential distortion  and  the  occurrence  of  rotating  stall  on  that  blade  row 
are  uncoupled  phenomena  except  under  conditions  where  the  disturbance  velocity 
approaches  the  propagation  speed  of  naturally  occurring  rotating  stall. 

This  matching  of  distortion  and  stall  propagations  speeds  is  not  likely  to 
be  attained  in  a typical  compressor  since  the  distortion  is  usually  stationary 
while  rotating  stall  propagates  at  some  fraction  of  the  rotor  speed.  The 
observed  dependence  of  stall  inception  on  the  extent  of  stationary  distortion 
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cells  (References  3 and  4)  would  appear  to  be  a result  of  phenomena  other 
than  the  time  to  stall  argument  (Reference  5).  The  dependence  may  be  due  to 
the  effects  of  a downstream  diffuser  such  as  studied  in  Reference  6. 


E.  CONCLUDING  REMARKS 

The  results  of  two  separate  experimental  investigations  of  rotating 
stall  have  been  presented  in  Section  II.  In  the  first  investigation, 

Section  II-B,  the  experiments  were  designed  to  investigate  whether  inter- 
ference between  a rotor  and  a stator  has  an  effect  on  the  conditions  for 
inception  of  rotating  stall.  The  rotor  row  and  the  stator  row  used  in  these 
experiments  had  previously  been  tested  in  isolation  so  that  the  conditions 
under  . iich  the  individual  blade  rows  undergo  inception  were  known.  Thus  the 
presence  of  interference  effects  between  the  two  rows  when  combined  to  form 
a compressor  stage  could  be  determined  by  comparison  with  the  individual  blade 
row  performances.  In  the  second  investigation,  Section  II-C,  the  effects  of 
circumferential  inlet  distortion  on  the  inception  and  properties  of  rotating 
stall  on  a rotor  were  studied.  These  tests  were  performed  both  with  the  dis- 
tortion pattern  stationary  in  the  laboratory  frame  of  reference  and  with  the 
distortion  pattern  rotating  about  the  compressor  axis  at  various  speeds. 

The  rotating  distortion  experiments  were  undertaken  to  determine  if  the  time 
of  passage  of  a rotor  blade  through  a distortion  cell  has  an  influence  on  the 
inception  of  rotating  stall.  As  with  the  blade  row  interference  studies,  the 
rotor  used  in  these  tests  had  been  tested  previously  in  isolation  so  that  the 
conditions  which  cause  rotating  stall  on  the  rotor  with  uniform  inflow  were 
available  for  comparison. 

In  all  of  the  rotor-stator  interference  experiments  of  Section  II-B, 
the  inlet  flow  conditions  relative  to  each  blade  row  were  such  that  in  the 
absence  of  interference,  rotating  stall  would  occur  first  on  the  rotor. 

These  studies  were  performed  at  a fixed  setting  of  the  rotor  blade  stagger 
angle,  but  with  four  different  settings  of  the  stator  blade  stagger  angle. 

The  combined  results  showed  an  increasingly  larger  delay  in  rotating  stall 


40 


inception  in  comparison  to  the  isolated  rotor  as  the  stator  stagger  angle  was 
increased,  that  is  as  the  stator  row  was  unloaded.  Thus  it  is  concluded  that 
rotor-stator  interference  does  have  an  effect  on  rotating  stall  inception  and 
that  the  magnitude  of  the  effect  depends  on  the  loading  of  the  blade  row  which 
would  be  unstalled  in  the  absence  of  interference.  The  results  of  these  ex- 
periments are  compared  with  the  predictions  of  two-blade  row  rotating  stall 
theory  in  Section  III. 

The  studies  of  rotating  stall  on  a rotor  in  the  presence  of  sta- 
tionary circumferential  inlet  distortion,  Section  II-C-2,  used  two  symmetric 
multi-lobe  distortion  screens.  Each  screen  covered  the  same  total  area  but 
one  consisted  of  two  lobes  and  one  consisted  of  four  lobes.  The  screens 
produced  square  wave  circumferential  inlet  distortion  patterns  with  a total 
pressure  drop  of  approximately  1.1  times  the  dynamic  pressure  far  upstream. 
Comparison  of  the  steady  state  turning  and  loss  performance  of  the  rotor 
with  and  without  inlet  distortion  showed  only  very  small  differences  in  these 
performances.  The  largest  difference  was  a small  constant  difference  in 
flow  turning  performance  between  the  isolated  rotor  and  the  rotor  with  inlet 
distortion.  Similar  comparisons  of  rotating  stall  inception  conditions  and 
ofthenumbers  of  cells  and  propagation  velocities  after  inception  showed 
negligible  differences  between  the  distorted  and  undistorted  flow  results. 
Rotating  stall  inception  and  its  properties  after  inception  appeared  to  be 
almost  completely  independent  of  the  presence  of  stationary  inlet  distortion. 

The  studies  of  rotating  stall  in  the  presence  of  rotating  inlet 
distortion,  Section  II-C-3,  were  performed  with  the  two-lobe  distortion 
screen  from  the  preceding  tests.  In  these  tests  the  distortion  screen  was 
rotated  about  the  compressor  axis  at  various  angular  velocities,  both  in  the 
direction  of  the  rotor  velocity  and  in  a direction  opposite  that  of  the  rotor 
velocity.  The  maximum  screen  rpm  in  both  directions  was  equal  to  the  mag- 
nitude of  the  rotor  rpm.  For  zero  and  negative  screen  rpm's  (screen  rotating 
in  a direction  opposite  that  of  the  rotor),  rotating  stall  inception  and  its 
properties  after  inception  appeared  to  be  almost  completely  independent  of 


the  presence  of  inlet  distortion.  However  the  resuHs  obtained  with  positive 
screen  rpm  were  greatly  dependent  on  screen  rotation.  When  the  screen  was 
rotated  at  a velocity  near  that  of  the  propagation  speed  of  rotating  stall 
on  the  isolated  rotor,  inception  occurred  considerably  earlier  and  rotating 
stall  appeared  to  occur  in  phase  with  the  distortion  screen  disturbances. 

At  screen  speeds  other  than  that  corresponding  to  the  natural  stall  propaga- 
tion speed,  rotating  stall  inception  was  delayed  relative  to  the  in-phase 
value.  When  the  compressor  flow  coefficient  at  inception  is  plotted  as  a 
function  of  screen/rotor  speed  ratio,  the  resulting  curve  is  similar  in  shape 
to  the  amplitude  versus  frequency  curve  of  a damped  harmonic  oscillator  with 
its  resonant  frequency  at  the  natural  propagation  speed  of  rotating  stall. 

The  overall  results  of  the  studies  with  stationary  and  rotating 
inlet  distortion  suggest  that  rotating  stall  inception  is  significantly  in- 
fluenced by  circumferential  inlet  distortion  only  when  the  distortion  dis- 
turbances are  nearly  in  phase  with  the  rotating  stall  which  would  occur 
naturally  on  the  same  rotor.  This  matching  of  distortion  and  rotating  stall 
is  not  likely  to  occur  in  a typical  compressor  with  stationary  inlet  dis- 
tortion. Experiments  which  display  a dependence  of  stall  inception  conditions 
on  the  circumferential  extent  of  stationary  distortion  require  an  explanation 
other  than  the  hypothesis  that  rotating  stall  depends  on  the  rotor  blade  pas- 
sage time  through  a distortion  cell.  Such  dependence  may  be  a direct  result 
of  the  geometry  of  the  distortion  pattern  in  combination  with  the  presence  of 
a diffuser  downstream  of  the  compressor  (Reference  6).  In  the  current  studies, 
diffuser  effects  were  minimized  by  extending  the  constant  area  annulus  far 
downstream  of  the  rotor  and  terminating  the  constant  area  annulus  with  an 
annular  honeycomb  structure. 

In  addition  to  the  results  discussed  above,  the  experimental  in- 
vestigations provided  qualitative  and  quantitative  data  which  are  valuable 
for  guidance  of  the  theoretical  development  and  for  assessing  the  success 
of  the  theoretical  predictions  by  comparison  with  experimental  results.  (See 
Section  III). 
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SECTION  III 

THEORETICAL  ROTATING  STALL  RESEARCH 

A.  INTRODUCTION 

The  primary  questions  that  have  been  addressed  by  the  theoretical 
studies  are  the  effects  of  distortion,  close  coupling  of  blade  rows  and 
compressibility  on  the  inception  conditions  for  rotating  stall.  The  previous 
program  at  Calspan  (reported  in  Reference  1)  showed  that  a two-dimensional 
small  disturbance  stability  theory  was  capable  of  adequately  predicting 
inception  conditions  for  isolated  blade  rows  of  high  hub-to-tip  ratio  in  low 
speed  flows.  Consequently,  the  major  portion  of  effort  during  the  present 
program  was  concerned  with  extending  this  theory  to  include  the  above  mentioned 
effects.  In  addition,  the  development  of  a three  dimensional  rotating  stall 
theory  was  initiated  and  its  present  state  of  development  is  reported  here. 

The  analysis  of  distortion  has  been  treated  by  an  extension  of  the 
small  disturbance  stability  theory  given  in  Reference  1.  This  strictly  limits 
considerations  to  small  distortions  of  the  mean  flow,  but  proves  adequate  to 
explain  the  qualitative  trends  found  in  the  experiments  reported  in  Section  II 
where  the  distortion  was  not  small.  The  distortion  analysis  is  presented  in 
Section  III-B. 

The  two-dimensional  stability  theory  had  been  formulated  for  a two 
blade  row  configuration  in  Reference  1.  A limited  amount  of  correlation  of  the 
two  blade  row  theory  was  also  given  there  but  tne  spacing  between  blade  rows 
for  the  configurations  analyzed  was  evidently  so  large  that  no  significant 
effects  of  spacing  on  rotating  stall  inception  were  found.  The  primary  effect 
found  for  those  configurations  was  that  the  blade  row  spacing  determined  the 
number  of  stall  cells  that  occurred  at  inception,  but  the  relative  swirl  angle 
for  inception  was  essentially  unaltered  from  the  isolated  blade  row  case.  The 
results  of  the  present  study  of  closely  coupled  blade  rows  drastically  alter 
this  second  conclusion.  The  present  results  indicate  that  for  a normal  type  of 
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closely  coupled  rotor-stator  stage,  the  inlet  swirl  to  the  rotor  at  which 
rotating  stall  occurs  is  considerably  higher  than  for  the  case  of  the  rotor 
alone.  The  analysis  for  blade-row  interference  is  presented  in  Section  III-C. 

The  previously  developed  stability  theory*  for  incompressible  flow 
has  also  been  extended  to  consider  the  case  for  wholly  subsonic  flow  external 
to  the  blade  row.  The  extension  proves  to  be  straightforward  and  does  not  turn 
up  any  new  mechanisms  that  influence  the  onset  of  rotating  stall.  Namely,  the 
slope  of  the  loss  curve  as  a function  of  inlet  swirl  is  still  the  dominant  blade 
row  characteristic  that  influences  the  stability  of  a given  flow  configuration. 
As  such,  the  theory  is  expected  to  give  an  adequate  prediction  of  the  rotating 
stall  boundary  provided  that  the  appropriate  loss  characteristics  of  the 
blade  row  in  compressible  flow  are  available.  This  analysis  is  presented  in 
Section  III-D. 


Although  a two  dimensional  theory  has  proved  adequate  to  describe  the 
rotating  stall  phenomena  that  have  been  observed  in  the  Calspan/Air  Force  .Annular 
Cascade  Facility  (a  low  speed  test  rig  of  high  hub-to-tip  ratio)  there  are 
man;.-  practical  configurations  where  this  situation  does  not  obtain.  In  particular, 
the  first  stages  of  many  actual  jet  engine  compressors  are  of  relatively  low 
hub-to-tip  ratio  where  three  dimensional  effects  may  be  important.  For  this 
reason,  the  development  of  a three  dimensional  stability  theory  was  initiated 
under  the  present  program.  The  overall  format  of  'he  theory  will  be  the  same 
as  in  the  two-dimensional  case  in  that  the  stability  of  a basic  mean  flow  to 
unsteady  perturbations  will  be  studied.  The  first  task  in  developing  such  a 
theory  is  to  obtain  a realistic  but  tractable  representation  for  the  mean 
steady  three  dimensional  flow  through  a blade  row  such  that  the  perturbation 
analysis  does  not  become  so  complex  as  to  obscure  the  physics  of  the  flow 
process.  A simplified  actuator-disc  model  has  been  developed  as  a viable 
candidate  for  this  required  mean  flow  representation.  The  model  is  not  limited 
to  small  turning  through  the  blade  row  and  uses  two-dimensional  blade  row  loss 
and  turning  performance  (such  as  would  be  available  from  conventional  cascade 
tests)  in  a strip  theory  fashion.  This  analysis  for  three  dimensional  flow  is 
presented  in  Section  II1-E. 
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B.  EFFECTS  OF  INLET  DISTORTION 

The  analysis  of  inlet  distortion  has  been  made  within  the  framework 
of  a two-diaensional  incompressible  flow  model.  It  was  shown  in  Reference  1 
that  this  type  of  analysis  was  adequate  for  high  hub-to-tip  ratio  configurations. 
Strictly  interpreted,  the  two-diaensional  assumption  restricts  the  analysis  to 
consideration  of  distortion  patterns  that  have  only  circumferential  variation. 
However,  it  was  found  in  Reference  1 that  some  radial  variations  in  flow 
quantities  can  be  handled  by  using  appropriate  radial  averages.  Presumably, 
this  would  also  be  true  for  the  distortion  problem. 

The  flow  aodel  employed  is  similar  to  that  used  in  the  small  disturbance 
stability  theory  of  Reference  1.  We  briefly  restate  the  basic  assumptions  that 
are  inherent  in  this  type  of  modeling: 

1.  The  flow  is  incompressible  (low  speed) 

2.  The  blade  row  has  high  hub-to-tip  ratio 

3.  The  annular  area  is  constant 

4.  The  flow  disturbances  have  wave  lengths  in  the  circumferential 
direction  ( Y -direction)  which  are  large  compared  to  either  the 
blade  spacing  or  chord. 

5.  The  local  performance  of  the  blade  row  is  only  a function  of  the 
local  inlet  swirl  and  dynamic  pressure  and  may  be  obtained  in  a 
quasi-steady  fashion  from  the  steady  state  performance  data. 

Assumption  4 in  conjunction  with  5 implies  that  the  admissible  dis- 
turbances have  a low  reduced  frequency  based  upon  thq  velocity  relative  to  the 
blade  and  the  blade  chord  length.  This  is  consistent  with  our  previous  results. 

It  has  been  suggested^  that  the  phenomenon  of  "stall  delay"  is  im- 
portant to  the  fashion  in  which  inlet  distortion  influences  the  stall  boundary 
of  a compressor  stage.  In  reality,  this  is  the  same  phenomenon  that  was  referred 
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to  as  a boundary  layer  lag  time  for  separation  in  the  early  rotating  stall 
works.  The  basic  idea  is  that  when  an  airfoil  (or  blade  row)  has  its  angle 
of  attack  (inlet  swirl)  suddenly  increased  to  a condition  where  the  airfoil 
would  normally  be  stalled  then  the  boundary  layer  does  not  respond  immediately 
to  the  new  conditions  but  requires  some  time  to  adjust.  During  this  adjustment 
period,  the  airfoil  (blade  row)  may  exceed  its  normal  steady  state  performance 
limits.  It  was  found  in  Reference  1 that  the  occurrence  of  this  phenomenon 
was  not  essential  to  the  development  of  a rotating  3 t a 1 1 ; consequently,  it  has 
not  been  considered  in  the  ensuing  analysis  of  distortion. 

In  addition  to  the  previously  mentioned  restrictions  on  the 
flow  model,  the  following  analysis  will  be  limited  to  the  case  where  the 
distortion  is  a small  perturbation  on  the  mean  steady  flow.  It  was  found 
necessary  to  limit  the  theory  in  this  respect  in  order  to  obtain  analytical 
results.  Moreover,  the  use  of  this  assumption  and  the  theory  of  linear  systems 
would  allow  us  to  immediately  state  the  results;  however,  since  the  linearised 
theory  adequately  explains  many  of  the  experimental  observations  given  in 
Section  1 1 -C , the  formal  development  will  be  presented  in  detail. 

The  flow  through  an  isolated  blade  row  in  incompressible  flow  is 
considered.  A finite  thickness  two-dimensional  actuator  as  shown  in  Figure  29 
is  used  to  model  the  blade  row.  Also,  an  arbitrary,  but  prescribed,  distortion 
in  total  pressure  and/or  velocity  is  indicated  ahead  of  the  blade  row . Sub- 
scripts o's  indicate  quantities  in  a blade  fixed  coordinate  system.  The  trans- 
formation between  blade  fixed  coordinates  and  the  laboratory  fixed  system  is 
given  by: 

x.  - x 

’ i - W6t 

t.  - t 

where  Wb  is  the  blade  velocity  in  the  laboratory  fixed  system.  The  conditions 
upstream  of  the  blade  row  are  denoted  by  a subscript  1 while  conditions  downstream 
of  the  blade  row  arc  denoted  by  a subscript  2.  The  mean  flows  upstream  and 
downstream  of  the  blade  row  are  uniform  but  of  different  swirl  angles.  The 
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mean  absolute  swirl  in  each  flow  region  is  denoted  by  <5,,  and  is  equal  to 
wy  (Jt  , The  blade  row  is  moving  with  an  absolute  .'elocity  and  the  non- 

dimensional  blade  velocity  is  given  by  /?  s */bA/(  . The  relative  swirl  with 
respect  to  the  blade  row  is  denoted  by  ^ and  is  given  by  r St  - h 
The  total  flow  quantities  which  are  denoted  by  a hat  are  decomposed  in  steady 
and  unsteady  parts  as: 

A 

VJ-  * W;  + ujt  (x  , y , t ) 

U i - Ul  * a-  ( x,  y , t ) 

- W,  t A t (X,  y.  t) 


where  H is  the  total  pressure. 


It  is  assumed  that  the  steady  s*ate  performance  of  the  blade  row 
is  known  in  the  following  form: 

The  turning  performance  of  the  Made  row  is  ;..en  by 
Jj  ■ G ( Jt ) 

and  the  losses  through  the  blade  row  are  giver,  ov 

r - / . * 


2 \ P ( Jo,  * *•,  I 

where  \ - \ f *f,  ) is  known. 


By  the  quasi-steady  assumption,  the  .nstantareows  values  of 

■i  H0.  in  unsteady  flow  are  obtained  *n  input  mg  the  required  ;r.star.t««necu» 

* *k  - - 

values  of  • and  W0i  in  the  above  functions.  *e  assume  that  the 

distortion  wave  does  not  alter  the  mean  part  of  these  quasi -steady  relations. 

In  the  actual  use  of  these  relations,  they  will  be  linear: :ed  about  their 
steady  state  values. 


We  start  by  considering  the  distortion  pattern.  Clearly,  the 
distortion  pattern  must  satisfy  the  equatio-s  of  motion.  Since  the  flow  is 
incompressible,  the  distortion  pattern  may  be  described  l a stream  function 


r 


ly(x,yt:)  which  is  related  to  the  velocity  components  by  U « 1y/iy  and 
W=-  Jfr/Jx  . The  equation  for  the  stream  function  is 


in 
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(1) 


where  H * - and  is  the  vorticity.  We  should  note  that  Equation  1 is  not 
limited  to  small  disturbances,  and  applies  to  a general  rotational  unsteady 
flow.  Now  we  are  interested  in  solutions  of  Equation  (1)  that  represent  wave 
like  patterns  superimposed  on  a uniformly  swirled  steady  mainstream.  That  is, 
we  seek  solutions  of  the  form 


1 f - U ■ w - W4  x ♦ a e 


j(O.X*by*Ct) 


j (I*  ♦ * n t ) 

ft  e (2) 


T st  two  terms  represent  the  mean  flow  with  swirl  » —■  and  the 

; ecw..«  two  terms  correspond  to  unsteady  waves  of  arbitrary  amplitude  OL  and 
6 respectively.  The  a.  b,  c,  1,  a and  n are  complex  constants.  Now  it  may 
be  verified  by  substitution  that  Equation  (1',  has  solutions  of  the  desired 
fora  provided  that  the  following  conditions  held: 
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Conditions  3 and  4 indicate  that  the  waves  are  converted  :r.  the 
mean  flow  direction.  Since  these  conditions  do  not  involve  the  wave  amplitude, 
Ol  and  .4  , we  arrive  at  the  very  important  conclusion  that  with  regard  to 
solution  forms  of  toe  type  C the  flow  acts  as  a linear  «v$te«.  It  is  not 
demonstrated  her*,  but  these  types  of  relations  general i:e  to  an  arbitrary 
number  of  waves  per  imposed  on  the  mean  flow.  Horeover.  *:  may  he  shown 
that  waves  which  are  different  narmomes  of  the  same  fundamental  frequence 
satisfy  Equation  ‘S'  . It  should  also  be  noted  that  the  disturbance  stream 
functions  employed  .n  Reference  1 are  of  the  form  considered  here. 
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These  results  have  the  following  significance  for  the  present  analysis 


1)  The  analysis  is  limited  to  small  distortion  because  the  blade 
row  performance  is  linearized  about  its  mean  value  not  because 
of  the  type  of  velocity  perturbations  considered. 

2)  The  usual  methods  of  Fourier  analysis  may  be  used  to  build  up  an 
arbitrary  distortion  wave  ahead  of  the  blade  row.  In  the 
following  work,  we  will  consider  a single  arbitrary  Fourier 
component  of  distortion  ahead  of  the  blade  row. 

3)  The  usual  method  of  studying  the  stability  of  a flow  is  to 

perturb  the  flow  with  a small  periodic  disturbance  and  determine 
under  what  conditions  the  disturbance  is  damped  or  amplified. 
When  the  basic  flow  is  time  dependent  as  in  the  present  case, 
this  type  of  analysis  can  become  quite  complicated.  However, 
since  the  perturbations  to  be  considered  are  of  the  same  form  as 
the  basic  Fourier  component  of  the  distortion  wave  ar.d  we  have 
just  shown  that  they  do  not  interact,  and  since  the  blade  row 
(actuator)  performance  is  linearized  its  mean  performance 

point,  we  may  formally  split  the  problem  into  two  parts.  The 
first  part  turns  out  to  be  the  stability  analysis  of  the  mean 
flow  without  the  distortion  wave:  This  is  precise!.-  th<  pror-.ee 
considered  in  Reference  i-  The  second  part  of  the  problem  :s 
the  response  of  the  blade  row  to  the  distortion  wave,  since 
the  first  part  ol  the  problem  is  adequately  treated  . n sef  ere-' . t 
we  will  concentrate  on  the  second  part  m the  following  wort 

Of  course,  the  same  split  up  of  tie  problem  can  be  accomplished 
by  a formal  linearization  procedure  whereby  we  would  consider 
the  stability  of  small  periodic  disturbances  superimposed  ipor 
a small  distortion  wave  However,  the  present  discoveries 
allow  this  split  up  without  recourse  to  arguments  that  . mo..t 
the  relative  orders  of  smallness  of  these  di sturbances . 


In  the  following  analysis,  we  will  consider  the  response  of  the  blade 
row  to  a distortion  wave.  It  suffices  to  consider  the  velocities  associated 
with  a distortion  wave  since  if  they  are  known,  the  definition  of  total 
pressure  aay  be  used  to  calculate  the  corresponding  total  pressure  variation 
of  the  wave.  The  distortion  problea  has  aany  similarities  to  the  stability 
problem  considered  in  Reference  1.  Namely,  the  solutions  for  the  disturbance 
velocities  upstream  and  downstreaa  of  the  blade  row  are  of  the  saae  general 
fora  and  contain  three  undetermined  constants.  These  constants  are  determined 
by  applying  the  saae  aatching  conditions  across  the  blade  row.  The  notation 
of  Reference  1 will  be  followed  as  closely  as  possible. 


It  is  sufficient  to  consider  the  response  of  the  blade  row  to  a 
single,  arbitrary,  Fourier  component  of  an  upstreaa  distortion  wave.  The 
perturbation  streaa  function  in  each  flow  region  is  readily  a.ailable  from 
the  work  of  Reference  1 as 
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where  the  , 8^  , are  constants  in  each  flow  region  4.  and  the  sub- 

script  n denotes  that  we  are  considering  the  n Fourier  component . Prescribing 
the  distortion  is  equivalent  to  specifying  D and  since  the  system  is  linear 
and  is  being  forced  at  a prescribed  frequency,  the  response  will  occur  only 
at  this  prescribed  frequenc  Hence,  we  need  only  consider  the  single  Fourier 
component  m each  flow  region  that  corresponds  to  this  forced  frequency. 

Since  the  Fourier  component  A the  distortion  is  specified,  then  j is  also 
specified  as  opposed  to  t>eir.g  a free  parameter  1..  the  stabilit>  analysis. 
Furthermore,  C will  be  all  re*.,  otherwise  the  spec i tied  distortion  wave 
would  grow  or  die  out  exponer*  all*  with  time.  Then,  the  overall  problem  will 
be  to  solve  for  the  response  . e . s , 3,^  *nd  34<<  for  a specified 
input.  3 . and  determine  conditions  for  the  response  to  become  unbounded 


The  remaining  bomwk ry  conditions  and  matching  conditions  will  be 
those  used  in  the  analysis  of  Reference  l.  In  the  following  analysis,  the 
subscript  for  Fourier  comp-nest  rs-jeber  .1.;  be  dropped  The  c rsditior.  of 
bounded  ness  at  upstream  ano  downstreaa  mfir.it-  requires  h .•  • 9 


>0 


Then 


assuming  that  the  distortion,  i.e.,  D,  , is  prescribed,  the  analysis  on  pp. 
46-53  of  Reference  1 can  be  extended  to  include  the  D,  terms.  Again,  three 
equations  are  required  to  determine  d,  , and  in  terms  of  D , . The 
same  three  matching  conditions;  namely,  mass  flow  conservation,  vorticity  com- 
patibility, and  the  flow  deflection  relation,  are  used.  With  these,  the 

homogeneous  system  of  equations  in  the  stability  theory,  Equation  (26)  of 
* 

Reference  1,  is  replaced  by  the  corresponding  inhomogeneous  system  of  equa- 
tions of  inlet  distortion  response  theory;  namely, 
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Equation  2*  of  Reference  ontam* 

'3,' 

in  Reference  1 should  read 
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Equation  (7)  is  a system  of  three  simultaneous  linear  equations  in 
B,/D,  , R2/D,  and  02/D,  . The  solutions  are 
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where  the  determinants  are 
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The  determinant  *17  can  be  expanded  ">ut  and  written  m the  loirs 
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where 


New  examining  these  results,  it  can  be  concluded  that  8,  , Dx  and  4^ 
will  be  finite  unless  od  • 0 . From  Equation  16  this  can  only  occur  if 
or  A • Aa  . The  first  case  is  impossible  since,  as  noted  previously,  C is 
completely  real  for  specified  inlet  distortion.  The  second  case  would  require 
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and  ff.P  (*,)  « A * — - (19b) 

But,  the  condition  described  by  Equation  (19a)  can  only  occur  if  the  inlet 
conditions  to  the  blade  row  are  such  that  the  blade  row  is  on  the  rotating 
stall  boundary  with  no  distortion  piesent  and  Equation  (19b)  requires  that 
the  prescribed  distortion  have  the  propagation  speed  for  the  self  excited 
rotating  stall  wave.  Hence,  we  find  that  the  distortion  will  have  no  effect 
on  the  flow  unless  the  undistorted  flow  conditions  would  produce  rotating  stall 
and  the  distortion  lias  a Fourier  component  corresponding  to  a natural  rotating 
stall  wave  under  those  conditions.  In  effect,  the  system  is  then  in  resonance, 
that  is.  we  are  then  forcing  a linear  system  at  its  natural  frequency.  Of 
course,  in  the  vicinity  of  resonance,  we  expect  the  velocity  perturbations  to 
become  so  large  that  the  linearization  is  no  longer  valid.  In  this  region, 
the  nonlinearit les  wilitake  over  to  limit  the  amplitudes  of  the  aot.on.  rhen, 
based  upon  the  linear  analysis,  we  expect  that  distortion  will  not  affect  tne 
rotating  stall  boundarv  unless  the  blade  row  is  near  stall  and  t.ne  dist»rti>n 
has  a Fourier  component  which  would  correspond  to  the  natural  rotating  stall 
mode  of  the  blade  row.  Comparison  with  the  results  of  the  e>;;-enaents  «ith  . 
rotating  distortion  screen  see  Figure  CS  s.-.ow  that  this  is  the  case 


It  snould  be  pointed  cut  that  a linearized  snail  perturbation  anal-sis 
would  result  in  a similar  conclusion  for  more  complex  flow  situations  such  as 
multiple  blade  rows  and  variable  annulus  area  Namely,  that  the  distortion 
would  not  affect  the  stall  line  unless  the  distortion  rut  the  s*ste*  it. 
resonance.  That  is.  once  the  linearized  stability  boundary  has  beer,  determined 
for  a giver,  flow  case  using  the  proper  blade  sow  performance  and  ros^ndar-  con- 
ditions for  that  case*,  the  addition  of  small  distortion  to  toe  flow  it.  t.ne 
fashion  done  here  will  only  result  in  unbounded  blade  row  response  i:  t ~.t 


S 5 


distortion  wave  has  a Fourier  component  corresponding  to  the  natural  rotating 
stall  mode  for  that  particular  case. 

It  should  be  reemphasized  that  the  foregoing  analysis  and  conclusions 
assumed  that  the  distortion  wave  did  not  influence  the  mean  steady  performance 
of  the  blade  row  and  the  blade  row  performance  was  linearized  about  its  mean 
steady  operating  point.  Some  effort  was  spent  upon  trying  to  remove  this 
assumption  from  the  analysis  and  thus  arrive  at  a nonlinear  stability  theory. 
Unfortunately,  this  type  of  analysis  cannot  proceed  very  far  without  resulting 
in  the  requirement  for  numerical  solutions.  The  analysis,  however,  did  reveal 
the  following  points  which  aid  in  the  understanding  of  how  distortion  affects 
the  stability  of  the  flow. 

1.  For  nonlinear  blade  row  performance  we  cannot  formally  split  the 
problem  into  two  parts,  i.e..  the  stability  and  response  problems 
cannot  be  treated  separately. 

2.  In  tne  nonlinear  case  the  response  of  the  actuator  contains  all 
the  Fourier  components  in  the  upstream  and  downstream  disturbance 
velocities.  This  results  in  the  distortion  wave  interacting  with 
the  blade  row  to  produce  a perturbation  of  the  aeui  steady  swirl 
upstream  and  downstream  of  the  blade  row.  This  perturbation 
dies  out  exponentially  with  distance  from  the  blade  row;  but  in 
effect,  results  m the  blade  row  seeing  a different  nean  inlet 
swirl  and  producing  different  mean  performance  cur  es  for  the 
blade  row  l.n  the  real  nonlinear  case  ther..  ea->.  individual 
distortion  wave  will  produce  new  mean  flow  conditions  for  the 
blade  row  which  arc  dependent  upon  the  type  of  distortion  being 
cons  t iered 

On  the  basis  of  the  theore:;%a«  analysis  performed,  we  were  unable  to 
draw  an*  general  conclusions  about  the  interaction  of  the  distortion  wave  and 
the  mean  flow,  but  on  the  basts  of  the  experimental  results  it  would  appear 
•hat  these  effects  are  oc.I»  important  near  resonance 
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THEORY  FOR  CLOSELY -COUPLES  BLADE  ROWS 


In  a typical  compressor  stage,  the  spacing  between  blade  rows  is  less 
than  a chord  length.  The  effect  of  this  close  spacing  has  been  investigated 
experimentally  in  the  annular  cascade  facility  using  a stage  consisting  of 
Rotor  Set  No.  1 followed  by  Stator  Set  No.  i.  The  experiments  are  described 
above  in  Section  II.  Here,  we  will  discuss  the  corresponding  theoretical 
calculations. 

The  two  blade  row  stability  theory  was  developed  and  applied  to  four 
cases  in  Reference  1.  These  cases  were  Stator  Sets  Nos.  1,  4 and  5,  each  m 
conjunction  with  the  stationary  inlet  guide  vanes,  and  Stator  Set  No.  6,  in 
conjunction  with  the  rotating  inlet  guide  vanes.  The  spacing  between  blade 
rows  was  large  (several  chord  lengths i in  all  of  these  cases  to  corresponds  to 
the  experimental  arrangement.  The  predicted  results  generally  agreed  well  with 
the  experiments.  One  case  was  rur  with  the  spacing  halved  and  the  results 
were  presented  in  Figure  SI  of  Reference  1.  Unfortunately,  there  were  no  data 
for  comparison 


In  the  present  study,  one  of  the  computer  programs  developed  tu  im- 
plement the  two  blade  row  analysis  of  Reference  l na<  Seen  modified  to  nondle 
the  input  tn  the  fora  m wnich  it  • « for  tne  st  »ge  cor. <ist  mg  >t 

Rotor  Set  No.  1 and  Stator  Set  V>  1 The  turning  and  loss  Jnr.i  for  Rotor 
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Jnfortanate  I • . tne  turning  and  io>s  dat  i tmt  *re  available  fir 
Stator  Set  No.  1 operating  m isolation  f igures  ,1  and  id  jf  Reference  I;  are 
limited  in  relative  inlet  swirl  angles  to  tne  range  of  angles  near  wfuc.n 
rotating  stall  xcjrre-J  on  this  set  The  r»  it . ■ - mle'  angles 
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Stator  Set  No.  1 in  the  present  experiments  are  those  corresponding  to  the 
outlet  swirl  of  Rotor  Set  No.  1 and  are  considerably  smaller  than  those  for 
which  Stator  Set  No.  1 data  exist.  Therefore,  extrapolations  of  the  data  were 
necessary.  The  way  in  which  the  data  have  been  extrapolated  is  described  in 
the  next  paragraphs. 

The  Stator  Set  No.  1 flow  turning  data  in  Figure  11  of  Reference  1 

m 

were  obtained  for  inlet  swirl  angles,  /S,  , greater  than  S3  degrees  because 
it  was  in  this  range  that  rotating  stall  occurred  on  this  blade  row.  For  the 

m 

present  calculations  these  data  must  be  extrapolated  down  to  /3,  values  of 
about  20  degrees.  The  extrapolation  was  done  as  follows.  First,  the  flow 
turning  data  for  Rotor  Set  No.  1 were  examined  because  they  had  been  obtained 
over  a wide  range  of  operation  (Figure  36  of  Reference  1).  Those  data  were 
recast  in  the  form  of  total  flow  turning  angle  |/3I(J  - versus  mean 

angle  of  attack  where  /3  and  are  the  inlet  and  outlet 

flow  angles  relative  to  the  blades,  respectively,  and  SgH  is  the  rotor  stagger 
angle.  The  results  are  plotted  in  Figure  31  (without  denoting  the  individual 
data  points),  for  <5,„  values  of  30,  40  and  SO  degrees.  These  data  collapse 
quite  well  into  a single  curve,  particularly  at  angles  of  attack  below  10 
degrees.  The  main  discrepancies  are  at  the  higher  angles  of  attack  for  each  <5»„ . 

Next,  the  Stator  Set  No.  1 flow  turning  data  (Figure  11  of  Reference  1) 
were  recast  in  the  same  form  of  | /3,  | - | J j/3,|  - j 6^  | , where  is 

the  stator  stagger  angle,  and  are  plotted  in  Figure  32  for  * 28.2,  3'. 2, 

and  S'. 2 degrees.  A single  correlation  curve  has  been  drawn  on  Figure  32  to 
pass  approximately  through  the  data  points  and  to  have  the  sane  behavior  at 
the  lower  angles  of  attack  that  the  Rotor  Set  No.  1 data  had  in  Figure  31. 

This  correlation  curve  in  Figure  32  was  then  used  to  generate  the  desired  ex- 
trapolations for  each  of  the  values.  The  operating  range  of  Stator  Set 

No.  1,  when  used  in  conjunction  with  Rotor  Set  No.  1,  includes  conditions  tor 
which  /3,  is  less  than  <5,„  . For  these  conditions,  the  flow  turning  is 
assumed  to  be  symmetrical  about  /4,  * - The  final  extrapolations  so 

generated  are  shown  m Figure  33,  where  the  actual  data  prints  are  indicated 
by  the  symbols. 
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The  extrapolation  method  for  the  loss  data  could  not  be  established 
on  such  a firn  basis.  No  loss  correlations  analogous  to  those  of  Figure  31 
could  be  found  for  the  data  of  Rotor  Set  No.  1 (Figure  30).  Instead,  those 
data  could  be  used  only  as  a guide  for  the  extrapolations  for  Stator  Set  No.  1, 
which  are  presented  in  Figure  34.  In  particular,  the  minimum  value  of  the 

m 

total  pressure  loss  coefficient,  AC pr  , was  taken  as  0.2  for  dCH  ■ 28.2  and 

37.2  degrees,  as  was  Measured  for  Rotor  Set  No.  1 at  the  two  lower  stagger 

» 

angles.  A larger  value  of  the  minimum  was  used  for  • 57.2  degrees, 

analogous  to  the  rotor  result  for  ■ SO  degrees.  The  loss  data  also  were 
assumed  to  be  symmetrical  about  fi,  ■ , just  as  for  the  flow  turning 

extrapolations  of  Figure  33.  It  is  unfortunate  that  a better  loss  data  extra- 
polation could  not  be  found  because,  as  discussed  in  Reference  1,  the  slope 
of  the  loss  curve  is  the  key  parameter  in  the  stability  analysis. 

Nuaerical  calculations  based  on  the  two  blade  row  stability  theory 
were  carried  out  for  the  stage  consisting  of  Rotor  Set  No.  1 followed  by 
Stator  Set  No.  1 using  the  extrapolations.  The  rotor  stagger  angle  was 
* 40  degrees  to  correspond  to  the  experiaent.  The  computed  damping 
factor  for  the  stage  is  plotted  as  a function  of  the  relative  inlet  swirl  to 
the  rotor,  , in  Figure  3Sa  both  for  the  rotor  alone  and  for  the  rotor  with 
the  stator  at  ■ 28.2  degrees.  Rotating  stall  inception  occurs  when  the 

damping  factor  first  goes  to  :ero  and  the  experimental  inception  points  are 
shown  for  comparison.  Rotating  stall  occurred  with  one  stall  cell  both  ex- 
perimentally and  theoretically.  The  presence  of  the  stator  delays  stall  both 
in  theory  and  experiaent  by  roughly  the  same  amount. 

The  computed  propagation  velocity  of  the  one-celi  rotating  stall 
pattern  is  plotted,  in  the  absolute  coordinate  system,  in  Figure  35b  as  a 
function  of  jb,  . again  for  the  rotor  alone  and  with  the  stator  at  <$%„  ■ 28.2 
degrees.  The  computed  curve  is  dashed  beyond  the  theoretical  inception  point 
to  indicate  that  the  theory  is  not  strictly  applicable  beyond  initial  incep- 
tion. The  experimental  propagation  velocities  at  inception  are  shown  as  well 
in  Figure  3Sb. 


Numerical  calculations  were  also  carried  out  at  stator  stagger  angles 
of  37.2  and  57.2  degrees.  At  * 37.2  degrees,  there  was  a slight  addi- 
tional delay  in  the  predicted  rotating  stall  inception  point,  but  not  as  much 
as  was  found  experimentally.  The  predicted  absolute  propagation  velocity  de- 
creased still  further  which  is  counter  to  the  experimental  results.  However, 
the  loss  data  at  <S,H  * 37.2  degrees  is  much  more  limited  than  at  6SM  * 28.2 
degrees,  see  Figure  34,  so  that  there  is  not  as  much  confidence  in  the  extra- 
polation. At  ■ 57.2  degrees,  there  is  only  a small  range  over  which  the 
extrapolated  stator  data  overlaps  with  the  rotor  data  of  Figure  30.  These 
rotor  data  were  obtained  at  rotor  speeds  up  to  1200  rpm.  Since  rotating  stall 
inception  occurred  experimentally  on  the  isolated  rotor  at  1000  rpm,  all  the 
performance  data  between  1000  and  1200  rpm  were  measured  with  rotating  stall 
present  and  so  include  very  high  losses.  For  numerical  calculations  at 

* 57.2  degrees  with  the  two  blade  row  stability  theory,  the  rotor  per- 
formance data  in  rotating  stall  were  used  between  1000  and  1200  rpm.  The 
discontinuity  in  a loss  curve  at  inception  was  faired  out  tor  this  purpose. 
Rotating  stall  inception  on  the  stage  with  * 57.2  degrees  was  not  pre- 

dicted up  to  the  1200  rpm  limit.  This  agrees  with  the  experiments  where  in- 
ception did  not  occur  until  a rotor  speed  of  1275  rpm  was  reached. 

The  results  for  a closely  coupled  rotor-stator  row  are  encouraging 
since  the  basic  trends  of  a delay  in  rotating  stall  inception  and  a reduction 
in  the  stall  propagation  velocity  are  predicted  qualitatively.  However,  the 
large  extrapolations  which  were  necessary  in  the  steady-st3te  loss  and  turning 
data,  make  quantitative  comparisons  of  limited  value.  More  significant  cor- 
relations must  await  the  acquisition  of  steady-state  loss  and  turning  data  on 
the  isolated  stator  row  at  the  required  inlet  conditions  as  well  as  the 
measurement  of  the  performance  of  the  rotor  row  and  stator  row  when  combined 
in  a closel>  coupled  stage. 

D.  EFFECTS  OF  COMPRESSIBILITY 

The  flow  model  and  overall  method  of  analysis  employed  m the  stability 
theory  treatment  of  coapressible  flow  through  a blade  row  are  quite  similar  to 
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the  incompressible  analysis.  The  major  differences  are  the  introduction  of 
one  more  dependent  variable,  the  density,  and  a corresponding  additional 
matching  condition  across  the  actuator  which  is  derived  from  the  energy 
equation.  The  flow  model  is  shown  in  Figure  36  where  again  the  blade  row  is 
modeled  by  a finite  thickness  actuator.  The  flow  region  upstream  of  the 
blade  row  is  noted  as  region  1 and  that  downstream  as  region  2.  The  mean 
swirls  upstream  and  downstream  are  constant  but  of  different  values.  Quantities 
in  a blade  fixed  coordinate  system  are  denoted  with  a subscript  o.  The  trans- 
formation between  blade  fixed  and  duct  fixed  corrdinates  is  given  by 

X.  - X 

i*  • i - wb  * 

t.  - t 

where  Wb  is  the  blade  velocity  in  the  duct  fixed  system.  This  velocity  is 

* W k 

given  in  non-dimensional  form  by  ■ — where  i denotes  the  flow  region  and 
UL  is  the  mean  axial  flow  in  region  i . The  mean  absolute  swirl  angle  in 
each  region  is  denoted  as  A and  the  tangent  of  6 is  defined  as  the 
absolute  swirl.  The  relative  swirl  with  respect  to  the  blade  row  is  given  by 

A ’ Si 


Following  References  1 and  2,  all  of  the  total  flow  quantities, 
are  denoted  by  a hat,  are  decomposed  into  steady  and  unsteady  parts  as 

U • J * a : r ' 
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/v  * 


P - 


C • w ♦ J f 

* 

where  P and  O are  the  pressure  and  density,  respect :•« eiy  . and  the  . den  tes 
the  flow  region.  The  assumptions  on  the  unsteady  flow  are  oasicall.  tne  s use 
as  in  the  incompressible  case.  ie..  the  flow  between  Made  rows  is  cons.Jereu 
mviscid,  so  that  frictional  and  belt  conduction  effects  are  neglected. 
Generally,  however,  the  flow  will  -'e  rotational  to  reflect  tne  losses  in- 
curred m the  blade  row.  The  unsteady  disturbance  quantities  are  assumed  to 
be  much  smaller  than  the  steady  quantities  so  that  t.ne  equations  of  notion  »an 
be  linearized  about  the  steady  flow. 


VI 


The  resulting  linearized  equations  of  motion  are  those  of  continuity, 
x-  and  y-  momentum  and  energy;  namely. 
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where  the  speed  of  sound  a.,  is  given  by  a*  * and  is  constant  in  each 

flow  region;  / is  the  ratio  of  specific  heats.  It  will  be  observed  that  p, 
and  0 are  coupled  with  the  velocity  components.  This  is  in  contrast  to  the 
incompressible  case  where  pi  is  zero  and  decouples.  In  the  incompressible 
case,  a stream  function  can  be  defined  and  the  equation  solved  for  it,  where- 
upon j.  . jj  and  jot  can  be  evaluated.  A stream  function  does  not  exist  in 
unsteady  compressible  fl>w,  but  Equations  (20)  to  (23)  can  be  solved  by  assuming 
solutions  of  the  form 
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,<*rc  “ * -«*/  anJ  Is  number  of  stall  ceils  and  r is  the  mean 

radius  ot  the  blade  row.  The  resulting  disturbance  is  spatially  periodic  in 
the  ^ direction  with  period  Zrtn  r . The  flow  is  neutrally  stable  if  C4  - 0. 
unstable  if  " and  stable  if  C^’C. 
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where  /?;  , 0-  , 0±  , 5t  are  constants  to  be  determined  by  the  boundary  con- 
ditions and  the  matching  conditions,  the  axial^  Mt ch  number  is  Mi  ■ W* it 
* (c  * — \Z^)/ Ui  , and  ot/  [f  - M*  ( t + 
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Reference 


2!  *r  51  reduce  to  l '*■  i-csrifrtis^U  sslu*;?*  -i 
• if  v* . i*  se*  equal  tc  :ers  and  the  c.  21 sturt-srcr  is  neglected 
itokficr.  the  n wti»  • -*<  • s 4 * t f f r • fras  that  ^f  lefere^t  1-  Sc  ;r.  the  c»i- 
’resssbie  case.  us?  j*  n the  ;nc  vnsrres  * ; b . e „i»e.  there  are  :rrotit;w,»; 

*•.  snirbance*  ten**  proport:  n,  to  - and  5 and  vorticity  waves  con- 
victed along  the  near  t . jh  >:rfiH.,nes  tern  pr~  port  canal  to  ^ In 

aid  it ; on  . tor  the  compress inle  f 1 _>«  .isf , the  tera  proportional  to  £ repre- 
sents an  ent  ropv  nave  convev ted  along  th»*  scan  flow  streamline 


There  are  eight  constant*  to  he  oeteramej  Iron  the  soundarv  con- 
ditions The  boundarv  conditions  at  upstream  and  downstream  infinity  allow 
the  lnstediate  determination  of  four  of  these  constants  Requiring  irrotational 
isenstropiv  flow  apsrream  of  the  blade  row  gives  J:,  * f • O . Requiring 
hounded  solutions  at  upstream  and  downstream  infinity  gives  ",  1 3.  • 0 . The 
tour  remaining  constants  are  determined  hi  the  matching  conditions  across  the 
tetuator;  namely,  the  conservation  of  mas*  flow,  conservation  of  vorticity, 
flow  turning  relation  and  conservation  of  energy.  The  first  three  o*’  these 
matching  conditions  are  similar  to  the  incompressible  case,  with  considera- 
tion for  variable  density,  and  the  fourth  is  new  for  the  compressible  flow  case. 


These  matching  conditions  will  be  outlined  in  the  following.  The 
conservation  of  mass  flow  across  the  actuator  is 
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In  vector  form,  the  vorticity  compatabili ty  condition  is 
j—  / • d * - j),  (V  < fi)  • dr  « - ^ 7 p ) • d ■*  + 


d*  + & F -d*  (33) 


where  the  line  integrals  are  taken  around  the  countour  indicated  in  Figure  36 
and  F indicates  the  frictional  forces.  Introducing  the  total  pressure  loss 
coefficient  \ as  defind  in  Section  1 1 1 B and  taking  the  limit  as  side  2 
approaches  side  4 ‘iquation  33  becomes: 
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where  kt  is  the  coefficient 

of  thermal 

conductivity  in 

region  t . UR 

is  the 

gas  constant 

and  at  is  the  t 

ound  speed 

in  regisii  <.  . 

In  the  above  expressions 

the  integrals  appearing  are  calculated  assuming  that  the  functions  appearing  in 
the  integrand  have  a linear  variation  between  blade  row  inlet  and  outlet  con- 
ditions. The  turning  relation  is  of  exactly  the  same  form  as  in  the  incom- 
pressible case,  Equation  (25)  Reference  1.  Substituting  Equation  (28)  through 
(31)  into  those  matching  conditions  results  in  the  following  homogenious  set  of 
equations . 
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The  characteristic  equation  is  then 


(37) 


which  gives  the  allowable  values  of  A (or  C ) for  self  excited  disturbance 
waves . 


Equation  (37)  is  not  a simple  polynomial  as  was  the  case  for  incom- 
pressible flow,  hence  we  cannot  make  an  analytical  comparison  with  that  case. 
However,  Equation  (37)  has  been  solved  numerically  for  several  cases  and  the 
calculations  for  Stator  Set  No.  4 (from  Reference  1)  at  a stator  stagger 
angle  of  28.2  degrees  are  presented  in  Figure  37  as  being  representative  of 

) 

| the  cases  examined.  In  these  calculations,  the  experimental  low  speed  turning 

l and  loss  performance  (given  in  Figures  12  and  15  of  Reference  1)  have  been 

l used  as  the  required  inputs  to  the  theory.  This  data  was  used  for  two  reasons. 


First,  Appropriate  5; *4*  ro»  perfsm^e  data  <adtr  corner* ssible  ils*  z or-i : 1 ion* 
•it  not  Available  Sacoodi*  , ^s«  of  in*  Is*  speed  dais  allows  a direct  ; '« • 
partsor.  with  tit*  results  of  ?Se  incompressible  theoretical  results  and  there 
by  isolate  the  effects  of  compressible  flow  alone  on  the  staM.it>  boundary 
figure  J*  indicates  that  compressible  flow  outside  the  blade  row  has  a slight 
destabilizing  effect  or.  relating  stall  inception  but  results  m a sizeable 
increase  m the  stall  propagation  velocity  relati.e  to  the  blade  row.  The 
slight  effect  on  inception  conditions  would  indicate  that  compressibility  does 
not  significantly  alter  the  flow  mechanism  involved  m rotating  stall.  More- 
over  the  slope  of  the  loss  curve  with  inlet  swirl  still  appears  to  dominate 
the  stability  of  a given  flow  configuration. 

The  present  results  suggest  the  possibility  that  the  stability  of 
a compressible  flow  configuration  may  be  approximated  by  using  the  compressible 
flow  blade  row  performance  in  the  incompressible  theory.  Concrete  proof  of 
this  contention  would  require  using  compressible  flow  blade  row  performance  in 
ooth  theories  and  comparison  of  the  results. 

E.  THREE-DIMENSIONAL  THEORY 

Although  the  various  versions  of  the  two-dimensional  small  disturbance 
stability  theory  have  been  reasonably  successful  in  predicting  the  rotating 
stall  boundaries  for  the  configurations  tested  in  the  Calspan/Air  Force  Annular 
Cascade  facility,  there  are  many  practical  compressor  designs  of  low  hub-to- 
top  ratio  where  the  two-dimensional  assumptions  are  not  expected  to  be  valid. 

This  has  been  the  motivation  for  the  three-dimensional  theoretical  studies. 
However,  the  overall  approach  of  studying  the  stability  of  small  disturbances 
superimposed  on  a mean  steady  flow  that  was  used  in  the  two-dimensional  case 
appears  to  be  substantiated  by  the  correlations  with  experiment  and  will  serve 
as  the  basis  for  the  three-dimensional  analysis.  The  first  step  in  this  approach 
for  the  three-dimensional  case  is  the  development  or  selection  of  an  appropriate 
analysis  for  the  mean  steady  portion  of  the  flow.  Once  this  portion  of  the 
flow  model  is  available,  an  appropriate  stability  analysis  may  be  performed. 
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Several  ttttedi  sf  prriicU5|  tVf  stead*  ? 1 o*  throwgi  .»  blade  ro%  rniitbU 
jn  tM  literature  »fr»  surveyed  but  non*  «N'r*  (ouvd  utitfKUrr,  for  wtaut 
natmi , tad  tSerefore . *n  appro^-iate  approximate  wtVxi  «*as  dev  el  ped  v 
brief  of  the  regu irwwti  ot  the  stead*  f lo*  model  and  of  atolls  al- 

ready available  will  precede  the  description  of  the  approuaate  model 
developed  under  the  present  program. 

It  is  required  that  the  model  for  the  mean  floe  be  capable  of  repre- 
senting the  losses  and  large  turning  that  occur  across  a blade  row  *hef>  it  is 
near  the  rotating  stall  boundary.  This  essentially  Units  consideration  to 
sone  sort  of  actuator  disk  model  where  enpincal  turning  and  loss  data  -say  be 
introduced.  Although  three-diaensional  mviscid  calculations  through  a blade 

g ) 

row  using  finite  difference  schemes  have  been  developed  by  Rae  * , it  is  not. 
yet  possible  to  introduce  viscous  or  separation  losses  in  these  calculation 
procedures.  As  such,  the  flow  through  the  blade  passages  cannot  be  calculated 
in  detail  but  must  be  modeled  by  an  actuator.  The  overall  approach  the.,  dua- 
lists of  developing  axisymetric  solutions  for  the  flow  upstream  and  downstream 
of  the  blade  row  and  matching  them  with  th«>  appropriate  conditions  across  the 
actuator  much  in  the  same  fashion  that  was  employed  in  the  two-dimensional 
theory.1  Such  a model  was  developed  by  Marble10  but  unfortunately  this  model 
is  limited  to  small  turning  through  the  actuator.  Actuator  models  allowing 
for  large  turning  through  the  blade  row  have  been  developed  by  Wu11  and  Oates.1 
These  works,  however,  use  an  eigen  function  expansion  procedure  which  result  in 
analytical  expressions  for  the  steady  velocity  components  which  are  too  com- 
plicated to  allow  for  a stability  analysis  for  the  unsteady  flow.  In  addition, 
these  approaches,  as  currently  formulated,  are  limited  to  flows  in  constant 
area  annuli.  Eigen  function  expansion  solutions  were  also  given  for  unsteady 
perturbations  superimposed  on  free  vortex  and  wheel  swirled  flows  in  constant 
area  annuli  by  Kerrebrock. ^ This  work  only  considers  the  propagating  wave 
forms  that  occur  in  such  flows  and  does  not  include  the  decaying  modes  that 

would,  in  general,  be  required  to  represent  a blade  row.  Also,  this  work 

1 4 

suffers  the  limitation  to  constant  area  annuli.  Oates  and  Knight-  and  Oates 

IS 

and  Carey  have  developed  a fairly  general  calculation  procedure  using  the 
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of  fssui  ilMNfl ts  These  Mtroit  would  be  technically  suitable  tor 

the  pretest  purposes  except  that  the  ru»ericai  calculations  involved  were 
beroad  the  scope  of  the  pretest  program.  The  actual  method  developed  aa>  be 
considered  at  a simplified  version  of  the  finite  element  approach. 


General  Theory 


The  present  analysis  considers  the  steads  flow  of  an  incompressible 
fluid  through  a blade  row  m a constant  area  annulus.  The  methods  employed, 
however,  will  leave  open  the  options  for  extending  the  analysis  to  compressible 
flow  through  a variable  area  annulus.  The  actuator  will  be  considered  as 
infinitely  thin,  since,  the  extension  to  finite  thickness,  is  not  considered 
difficult.  The  basi*  idea  of  the  present  analysis  is  to  represent  the  flow 
quantities  in  the  regions  upstream  and  downstream  of  the  actuator,  by  poly- 
nomials in  the  radial  variable  with  coefficients  that  are  functions  of  th-. 
axial  variable.  These  coefficient;  will  then  N»  d»r*min#H  hy  sari  trying  the 
equations  of  motion  in  an  integrated  or  average  sense. 


The  basic  system  of  momentum  and  continuity  equations  for  steady 
axisymetric  flow  are 
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The  coordinate  systems  is  indicated  in  the  following  sketen 


The  velocities  in  the  t , 9 and  2 direction  are  denoted  as  V , w and  U 
respectively.  The  static  pressure  is  denoted  as  P , The  actuator  is  located 
at  2*0.  The  annulus  of  constant  area  has  an  inner  radius  denoted  by  rM  and 
an  outer  radius  denoted  by  rT  . The  boundary  conditions  are: 


e No  flow  through  the  annulus  w*n« 

V - O at  r - rf  , rH 

• Uniform  flow  at  upstream  infinity 

U-^U,  , V,  W — 0 ; 2 — » - 


e All  quantities  bounded  at  downstream  infinity 


# Flow  turning  and  loss  relationships  specified  through  the  actuator 
e The  U and  V velocity  components  are  continuous  through  the  actuator 


Instead  of  working  directly  with  the  velocity  components  and  pressure 
as  dependent  variables,  it  was  found  to  be  more  convenient  to  work  with  the 
following  defined  quantities: 
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* A stress  function  defined  such  that 
• - ~7 


— 


then  Equation  J8d  is  satisfied  identical! 


The  total  pressure  divided  by  dens  it? 


O * 

The  product  r w 


r • rt 


e The  £?  component  of  vorticity 

ii  . *1  . >!L 

• Si  Sr 

Substituting  these  definitions  into  Equation  (38)  and  using  the 
definition  of  vorticity,  the  governing  equations  for  these  quantities  become: 


Kr  * - 

— V 

r ** 

m — 

m0 

(41) 

i h K n 

dr  r * 

r \»s 

r * 

3 

9 r 

( rus ) - o 

(42) 

32  r 9 

rus 

r* 

d 

32 

(r  u>  ) • 0 

(43) 

Vr  JL  (ruj)  - If,  J-  (rus)  - 0 (44) 


Upstream  of  the  blade  row  we  must  have  dl  s 0,H  « constant  and 
W « constant  ■ 0,  and  the  only  unknown  flow  quantity  is  y which  is  determined 
by  the  homogeneous  version  of  Equation  (41).  Now  y is  to  be  approximated  by 
a polynomial  in  f and  it  has  been  found  that  a three  term  polynomial  is  the 
lowest  order  polynomial  that  may  be  used  without  resulting  in  a trivial  solu- 
tion. Hence,  for  the  initial  stages  of  the  analysis,  we  consider  only  a three 
term  polynomial  approximation  for  the  flow  quantities  with  the  understanding 
that  no  steps  in  the  analysis  will  be  used  that  would  prevent  generalization 
to  higher  order  polynomials. 
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He  fine  define  a perturbation  stream  function  by 
V * -J  V.  r*  ♦ 0 

Then  Me  find  that  $ satisfies  the  following  equation. 

* <*>"  ' ~ br  • 0 (45) 

The  appropriate  $ will  be  of  fora 

<p  - a0(M)  * at(z)r  ♦ eLi(z)rJ  (46) 


The  a*  are  to  be  deterained  by  satisfying  Equation  (45)  on  an 
integrated  basis.  Naawly,  we  require  that 


' 7 ) 


dr 


(47) 


for  J • o,  > , ... 

Sufficient  values  of  ^ are  used  to  provide  a number  of  equations  equal  to  the 
number  of  unknown  functions.  Using  the  wall  boundary  condition,  Equation  (39) 
we  find  that 


ao  ■ 


where 


a,  . K,aA  ♦ 

K,  ‘ rr 

K,  “ ♦ '>> 


and  -4,  and  ■ki  are  constants  to  be  determined.  This  effectively  leaves  only 
one  unknown  function,  a4  » so  that  only  the  j ■ 0 version  of  Equation  (47)  is 
needed.  This  becomes  an  ordinary  differential  equation  with  constant  coefficients 
for  viz. 


R a" 
i i 


(48) 
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where 


(r 


i 


k 


Kt  JLn  — 

rT 


and  che  prices  denote  differentiation  with  respect  to  l . The  general  solution 
of  Equation  (48)  is 


0 - tnt 

a - 4 -£  ♦ c e 

X *1  fl  ' 


* C2e 


(49) 


/ fl  \'h 

where  and  is  real.  Application  of  the  upstream  boundary  condition 

gives 

4 


« C. 


and  Cx  remains  unknown  at  this  stage.  Cx  is  renamed  6>a  for  convenience  of 
the  following  analysis.  To  this  stage  of  the  analysis,  we  have  then 


- 

ml 

K.dx  e 

(50a) 

- 

ml 

Ki  e 

(50b) 

- 

ml 

(50c) 

Now  proceeding  to  the  downstream  flow,  we  will  denote  <p  downstream 
of  the  actuator  by  and,  in  general,  all  downstream  quantities  will  be 

denoted  by  a star.  The  partial  differential  equation  for  <p*  is 


*rr  * *1.  - y *r  - “ r il* 

Using  the  polynomial  approximation  for  0*  we  let 

, * * * * 4 


(51) 


(52) 
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The  wall  boundary  conditions  again  give 


a. 


< a. 

• i 


so  that  <1*  is  the  only  remaining  unknown  function  and  is  found  by  satisfying 
the  integrated  version  of  Equation  (4?);  namely, 


(53) 


where  r r» 

I,  •-/  r SI*  d r 
’■''V 

Before  proceeding  to  solve  this  equation,  a better  representation  of  I,  is 
necessary.  The  first  step  in  obtaining  this  representation  is  accomplished  by 
multiplying  Equation  (42)  by  y*  and  subtracting  Equation  (43)  multiplied 
by  1/*  . The  result  is 


« LH* 

r ii 


IN 
3 


(54) 


or  H « H Ufr  ) and  consistent  with  the  polynomial  approximation  procedure 


H ' - ft  V * ♦ D, 

where  ft  and  D,  are  constants. 


(55) 


Also  from  Equation  (44)  we  ’nave 

ruj*  - F (y*)  - BV*  + Dx  (56) 

where  3 and  D2  are  constants. 

Substituting  Equations  (56)  and  (55)  into  either  Equations  (42)  or  (43), 
results  in 
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V*  + o3 
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(57) 
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(58 


Using  this  equation  to  evaluate  If,  gives 

I,  * 2 m2  - B2 a* 


where 


«•  ’ A {»  ^TT  - B ' T ^ *’>  - Wr-V 


r 


Then  using,  Equation  (58)  in  Equation  (53),  we  obtain 


3 * 


Aa  - M a#  - 4 — 2-  + 2 m G, 


(59 


where 


..  a i .a 

M * m - B 


The  general  solution  to  this  equation  is 


a. 


• b.  e 


-Mi 
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where  bx*  0 to  keep  0*  bounded  at  infinity  downstream.  The  continuity  of 
U and  1/  across  the  actuatcr  is  equivalent  tc  keeping  <f>  and  <p'  continuous 
across  the  actuator  or 


and 


4>  (o)  m <p  k(0) 
4>'Co)  - <j)*'(o) 
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which  give 


and  " 


*z  ' 0 

2 m*  Q, 

M*  (Mi-ryn) 

- 2 m*Qt 
M (M  + m) 


Thus  far,  we  have 


(63) 


a,  - 


-2rr)XQt 


z m(M  + m) 

V - y + (Km+K,r  + r* ) 

f**  y U/*  + (K0  + K,r  * rx)  a* 


(64) 

(65) 


so  that  the  steady  velocity  field  is  known  up  to  the  four  constants  R , 0f  , B 
and  Dz  (appearing  in  Q1  and  M ) which  are  associated  with  the  steady  turning 
and  loss  performance  of  the  blade  vow.  In  the  previously  mentioned  actuator 
disk  theories,*0-15  the  turning  ( W component  of  velocity)  and  the  total 
pressure  are  considered  as  specified  at  the  outlet  of  the  blade  row  (actuator 
sheet)  which  is  equivalent  to  specifying  these  four  constants.  In  fact,  for 
the  present  purposes,  they  could  be  chosen  so  as  to  agree  with  the  annular 
cascade  measurements.  However,  for  the  unsteady  perturbation  analysis  that 
will  employ  this  model  for  the  mean  flow,  the  blade  row  performance  as  a 
function  of  local  inlet  swirl  will  have  to  be  used.  Therefore,  we  will 
develop  the  steady  analysis  using  the  appropriate  form  of  matching  conditions 
derived  from  the  local  blade  row  performance  data.  Before  presenting  an 
example  calculation,  a few  comments  concerning  the  form  of  the  solution. 
Equations  (63)  through  (65)  are  in  order. 

There  is  no  tangential  velocity  perturbation  transmitted  upstream  of 
the  blade  row.  There  is,  however,  an  axial  velocity  and  a radial  velocity 
perturbation  upstream  of  the  blade  row.  This  is  to  be  contrasted  with  the 
two-dimensional  steady  case  where  no  disturbances  are  transmitted  upstream  of 
the  blade  row.  This  feature  of  the  three  dimensional  theory  is  expected  to  be 
very  significant  in  considering  blade  row  interactions  that  are  currently  not 
explainable  by  two-dimensional  theories.  It  is  noted  that  for  the  specific 
polynomial  form  used  here,  a single  linear  differential  equation  is  obtained 
for  the  coefficients.  If  a higher  order  polynomial  had  been  used,  the  result 
would  be  a system  of  differential  equations,  but  the  system  would  still  be 
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linear.  Also,  the  restriction  to  constant  annulus  area  results  in  a system 
of  equations  with  constant  coefficients.  Consideration  of  variable  area 
annuli  will  result  in  a linear  system  with  non-constant  coefficients. 

2.  Example  Calculations 

As  an  example  of  the  above  theory,  the  turning  and  loss  performance 
of  the  stationary  guide  vane  row  of  the  Calspan/Air  Force  Annular  Cascade  has 
been  estimated.  The  blade  row  turning  and  loss  performance  required  by  the 
theory  were  obtained  from  two-dimensional  cascade  data  given  in  Reference  16. 
The  data  are  applied  in  a strip-wise  fashion  at  each  radial  station  assuming 
that  the  outlet  conditions  to  the  blade  row  are  only  a function  of  the  local 
inlet  swirl  and  dynamic  pressure  at  the  given  zadial  location.  The  local 
inlet  conditions  are  not  those  corresponding  to  the  undisturbed  flow  far  ahead 
of  the  blado  row,  but  are  calculated  as  part  of  the  solution. 


The  turning  and  loss  relations  across  the  blade  row  will  be  given  for 
the  general  case  first  and  then  specialized  to  the  guide  vane  case.  These 
relations  are  applied  in  a coordinate  system  fixed  to  the  blades.  Denote 
conditions  immediately  upstream  of  the  blade  row  by  a subscript  1 and  those 
immediately  downstream  by  a subscript  2.  Then  at  any  radial  station,  the  total 
pressure  change  normalized  by  the  density  is 


h2-h, 


r u)  (Wz-  W,) 


L_ 

P 


(66) 


where  uj  is  the  blade  row  rotational  velocity  and  L represents  the  losses 
and  it  is  assumed  that  they  are  given  by 

L ■ J X P jV  + (K  - rw)l| 

The  total  pressure  loss  coefficient,  X » *s  defined  as  in  the  previous  section 
of  this  report  and  it  is  assumed  that  it  is  available  from  two-dimensional 
cascade  data  as  a function  of  solidity,  stagger  angle,  inlet  swirl  and  airfoil 
section  shape.  Since  these  parameters  may  vary  along  the  span  of  a blade,  X 
will  in  general  be  a function  of  t . Equation  (66)  is  merely  the  Euler  turbine 
equation  with  the  added  consideration  of  losses. 
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Since  the  upstream  total  pressure  is  constant  and  assumed  known,  the 
previous  results  for  H.,  and  W0  may  be  substituted  into  Equation  (66)  to  give 

A yr  * (0)  + D,  - A/,  - lu(BW  *(o ) t D2)  - --  (67) 

P 


where 
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(->  r?co>)x 


r 


(68) 


liquation  (67)  may  be  differentiated  with  respect  to  r to  give 

A Y*  (o)  = uj3V*(o)  - — C69) 

r dr 

which  really  represents  the  vorticity  compatabil ily  relation  for  three- 
dimensional  flow. 


Let  the  local  relative  swirl  ahead  of  and  behind  the  actuator  be 
denoted  by  so  that 


9 W'  ~ r uj 

\A  . — - 

It  is  then  assumed  that  the  turning  performance  of  the  blade  row  is  available 
from  two-dimensional  cascade  tests  so  that  £ is  known  by 


61  G i 3, ) 


(70) 


where  G is  also  a function  of  solidity,  stagger  angle  and  airfoil  shape, 
which  in  general  will  vary  along  the  blade  span.  Using  the  previous  forms  of 
the  solution  obtained,  Equation  (70)  may  be  written  as 

B¥*(o)  + Dx  « rzuj  + If*  (o)  R (71) 


differentiating  this  with  respect  to  r gives 

& Y*(o)  - zno  + lf*r(o)C  + yr*(o)Gr 


(72) 
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Now  integrating  Equations  (67),  (69),  (71)  and  (72)  between  CT  and  r*H  gives 
a system  of  four  algebraic  equations  for  the  four  unknown  constants 

The  solution  of  this  system  will  be  demonstrated  for  the  case  of 
the  fixed  guide  vane  row  used  in  the  Calspan/Air  Force  Annular  Cascade  Facility. 
This  blade  row  has  a hub  to  tip  ratio  of  0.8,  a constant  solidity  of  1.17 
along  the  span  and  untwisted  blades.  The  stagger  angle  is  constant  along  the 
span  and  is  variable  approximately  from  18°  to  52°.  The  guide  vanes  had  an 
NACA  63-  (24A^  K^)  1C  guide  vane  profile  which  is  10  percent  thick.  These 
parameters  do  not  correspond  to  those  tested  in  Reference  16;  the  solidity 
falls  within  the  test  range  but  all  the  tests  results  are  for  a 6 percent 
thick  section.  Corrections  could  have  been  applied  to  the  turning  performance, 
however,  there  was  no  reliable  method  to  correct  the  loss  data.  In  addition, 
preliminary  comparison  of  the  loss  data  presented  in  Reference  16  with  that 
measured  in  the  annular  cascade  facility  indicated  a large  descrepancy  in 
levels.  This  is  presumably  due  to  the  large  differences  in  Reynolds  number 
between  the  two  tests  as  well  as  the  thickness  differences.  The  Reynolds 
number  in  the  annular  cascade  facility  was  9.5  x 104  where  as  those  in 
Reference  16  were  approximately  2.5  x 10S.  It  was  realized  at  the  outset  that 
these  differences  in  loss  level  would  lead  to  poor  correlation  between  the 
present  theory  and  the  annular  cascade  data,  but  the  numerical  calculations 
were  continued  in  order  to  check  that  the  theory  was  completely  formulated 
and  would  give  the  right  trends.  The  experimental  loss  and  turning  data 
corresponding  to  a solidity  ratio  of  1.0  were  used  in  the  theory. 

For  this  particular  case  with  oj  - 0 and  W,  = 0 there  is  no 
relative  swirl  ahead  of  the  guide  vane  row.  In  addition,  the  constant  blade 
row  geometrical  properties  along  the  span  greatly  simplify  the  application 
of  the  matching  conditions  and  the  integrated  forms  of  Equations  (72)  and  (71) 
become  respectively: 
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8 = Jt.Gr  (75) 
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The  integration  o£  Equation  (69)  and  (67)  give  respectively 

R X (r»  + rr'>  * 


and 


U 


D,-H, 

Ua 


2 f 3 r * 
’ H r 


9 


- T * »*?,*) 


where 


A - J- 


(rT  * rH)z  (rT3- rH’) 


iz 


<rr-rn) 


Jn  — - Z 


r - r~ 
t H 


ru 


R S — ~ + 2 K t~+4  + 

2 \rT  r«  rH 


(rH  i rTy  (rH  - rT)x' 


rH*rr* 


(74) 


(7S) 


(76) 


Substitution  of  Equations  (73)  - (76)  into  the  definition  of  give  a single 
algebraic  equation  for  CJ,  namely, 
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and  = " = .8  for  this  case.  Equation  (77)  has  been  solved  numerically 
and  has  been  found  to  possess  only  one  real  root-  The  required  G and  \ 
functions  at  each  guide  vane  stagger  angle  were  obtained  from  Reference  16  as 
previously  explained.  In  terms  of  these  parameters,  the  overall  average  total 
pressure  loss  coefficient  for  the  blade  row  is 


2 (D,  - Ht ) 


u; 


(78) 


This  expression  is  plotted  as  a function  of  stagger  angle  in  Figure  38  where 
the  corresponding  experimental  measurements  are  also  shown.  As  previously 
explained  because  of  the  difference  in  loss  levels,  poor  quantitative  agreement 
was  expected.  However,  throughout  most  of  the  stagger  angle  range  the  theory 
has  the  proper  trend.  Comparison  of  the  theoretical  results  with  the  two- 
dimensional  cascade  data  of  Reference  16  also  shown  in  Figure  38,  indicates 
that  the  theory  is  predicting  the  induced  drag  on  the  blade  row  correctly. 

The  reason  for  the  abnormally  high  loss  prediction  at  a stator  stagger  angle 
of  48°  is  not  clear.  Perhaps  a three  term  polynomial  is  inadequate  to  represent 
the  large  radial  variations  that  occui  under  these  high  loss  conditions. 
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The  theoretical  prediction  of  turning  performance  is  compared  with 
the  experiment  in  Figure  39.  The  theory  is  about  four  degrees  low  throughout 
the  stagger  angle  range.  This  is  easily  explained  by  the  difference  in  solidity 
and  thickness  between  the  annular  cascade  and  the  data  used  from  Reference  16. 
Overall  then,  it  is  seen  that  the  theory  predicts  the  qualitative  trends 
correctly.  The  absolute  level  of  the  predictions  is  off  because  the  two- 
dimensional  cascade  data  was  not  strictly  appropriate  for  the  annular 
cascade  guide  vane  geometry.  Moreover,  better  agreement  of  the  levels  of 
turning  and  loss  performance  between  theory  and  experiment  may  be  forced  by 
empirical  specification  of  the  four  unknown  constants  R,B,  D1  and  0 . 

Therefore,  the  present  theory  should  serve  as  an  adequate  representation  of  the 
mean  flow  upon  which  we  can  perform  an  unsteady  stability  analysis. 

F.  CONCLUDING  REMARKS  ON  THEORETICAL  INVESTIGATIONS 

During  this  program,  the  two-dimensional  stability  theory  was  ex- 
tended to  consider  inlet  distortion  for  incompressible  flows  and  compressible 
flows  without  inlet  distortion.  The  analyses  show  that  the  principal  effects 
of  distortion  occur  when  the  distortion  wave  corresponds  to  the  natural 
rotating  stall  mode  for  the  system.  The  two-dimensional  theory  was  also  used 
to  examine  the  stability  of  the  case  of  two  closely  coupled  blade  rows  in  a 
rotor-stator  stage.  The  theory  indicates  that  blade  row  interference  fc~  this 
configuration  has  a large  stabilizing  effect  on  the  stall  boundary-  The  cor- 
relations with  experiments  verify  this  trend  but  are  hampered  because  of  lack 
of  appropriate  steady  blade  row  performance  data. 

In  addition  to  these  two-dimensional  analysis,  the  initial  develop- 
ment of  a three-dimensional  theory  is  reported.  An  approximate  three- 
dimensional  analysis  for  the  steady  flow  through  a blade  row  is  given.  This 
model  is  to  serve  as  the  basis  of  a stability  analysis  to  predict  the  rotating 
stall  boundary. 


82 


SECTION  IV 

SUMMARY  AND  CONCLUSIONS 


A combined  experimental  and  theoretical  research  program  on  rotating 
stall  in  axial  flow  compressors  has  been  conducted.  The  primary  topics 
addressed  in  this  report  are  the  effects  of  stationary  and  moving  distortion 
on  the  rotating  stall  boundary  for  an  isolated  rotor  row  and  the  effects  on 
the  boundary  of  close  coupling  of  the  blade  rows  in  a rotor-stator  stage. 

The  tests  on  stationary  distortion  effects  on  a rotor  row  included 
both  a two  lobed  and  a four  lobed  distortion  pattern.  Although  the  total 
pressure  distortions  were  quite  high  for  both  patterns,  they  did  not  signifi- 
cantly affect  the  rotating  stall  boundary  or  the  rotating  stall  properties 
after  inception.  Similarly,  measurements  of  averaged  blade  row  loss  and 
turning  data  showed  little  influence  of  the  distortion.  The  tests  with  moving 
distortion  were  performed  by  rotating  the  two  lobed  pattern  about  the  com- 
pressor axis  both  in  the  direction  of  blade  row  rotation  and  opposite  to  it. 

It  was  generally  found  that  the  rotating  distortion  pattern  either  had  little 
effect  on  or  else  destructively  interfered  with  rotating  stall  formation 
except  when  the  pattern  was  rotated  at  speeds  near  the  natural  stall  propagation 
velocity  for  undistorted  flow.  Under  these  latter  conditions,  the  rotor  in- 
curred rotating  stall  at  much  higher  flow  coefficients  than  for  the  undistorted 
flow  case.  These  results  contradict  the  "time  to  stall"  arguement  that  has 
been  used  in  several  instances  to  explain  the  effects  of  multilobe  distortion 
patterns  on  compressor  stall  boundaries.  The  theoretical  analysis  are  in 
general  agreement  with  the  experimental  findings.  The  theory  concludes  that 
distortion  will  not  influence  the  rotating  stall  boundary  unless  the  distortion 
pattern  contains  a Fourier  component  which  corresponds  to  the  natural  rotating 
stall  mode  for  the  system. 

The  experiments  on  the  effect  of  close  coupling  on  the  rotating  stall 
boundary  for  a rotor-stator  stage  show  that  the  blade  row  interference  results 
in  a significant  delay  in  rotating  stall  inception.  The  rotor-stator  stage 
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stalls  at  a much  larger  inlet  swirl  than  the  isolated  rotor.  The  theoretical 
results  also  show  this  delay,  but  the  numerical  predictions  of  the  theory  are 
generally  inaccurate  due  to  lack  of  the  steady  state  performance  data  for  the 
isolated  blade  rows  in  the  swirl  angle  ranges  needed.  The  theory  requires  these 
steady  state  data  as  inputs. 

In  addition  to  the  above  mentioned  theoretical  studies,  which  were 
performed  with  an  incompressible  two-dimensional  small  disturbance  stability 
theory,  studies  were  performed  to  extend  the  two-dimensional  theory  to  include 
compressibility  and  to  initiate  the  development  of  a three-dimensional  theory. 
The  theory  for  compressible  flow,  which  assumes  that  the  flow  is  subsonic 
everywhere,  indicates  that  the  slope  of  the  loss  curve  as  a function  inlet 
swirl  still  controls  the  stability  of  a given  flow  configuration.  The  three- 
dimensional  theory  has  not  been  completed  at  this  time,  but  the  development 
to  date  is  given  in  this  report. 
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Figure  3 


DETAILS  OF  ORIGINAL  ROTATING  HUB  CONFIG 
URATION  IN  ANNULAR  CASCADE  FACILITY 
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Figure  4 

ROTATING  HUB  CONFIGURATION  MODIFIED 
TO  PROVIDE  INDEPENDENT  ROTATION  OF 
TWO  ROTOR  ROWS 
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Figure  5 NOTATION  FOR  ANNULAR  CASCADE  WITH  ROTOR-STATOR  STAGE 


J79  CASING  LENGTH,  26.88 
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Figure  7 STATIC  PRESSURE  RISE  ACROSS  ROTOR  SET  NUMBER  1 
ROTOR  STAGGER  ANGLE,  6 RM  = 40° 
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Figure  8 (Cont.)  STATIC  PRESSURE  RISE  ACROSS  ROTOR-STATOR  SET  NUMBER  1 
ROTOR  STAGGER  ANGLE,  5RM  = 40° 


SYSTEM  RELATIVE  TO  ROTOR 


J79  CASING  LENGTH,  26.98 


F»gure  1 1 ANNULAR  CASCADE  CONFIGURATION  USED  FOR  ROTATING  STALL  STUDIES  WITH 
STATIONARY  AND  ROTATING  INLET  DISTORTION 
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TOTAL  PRESSURE  COEFFICIENTS,  Cp  (r,  9), 
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Figure  16  CIRCUMFERENTIAL  DISTRIBUTIONS  OF  SWIRL  ANGLE,  (3-,,  BETWEEN  DISTORTION 
SCREEN  AND  ROTOR,  SCREEN  CONFIGURATION  B 
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Figure  17  (Cont.) 

TOTAL  PRESSURE  COEFFICIENTS  DOWNSTREAM 
OF  ROTOR,  DISTORTION  SCREEN 
CONFIGURATION  A 
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Figure  18  (Cont.) 

TOTAL  PRESSURE  COEFFICIENTS  DOWNSTREAM 
OF  ROTOR, DISTORTION  SCREEN 
CONFIGURATION  B 
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Figure  19  CIRCUMFERENTIAL  DISTRIBUTION  OF  SWIRL  ANGLE  J DOWNSTREAM  OF 
ROTOR  SCREEN  CONFIGURATION  A 
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Figure  20  EFFECT  OF  DISTORTION  SCREEN  CONFIGURATION  A ON  AVERAGE  SWIRL  ANGLES 
DOWNSTREAM  OF  ROTOR  SET  NO.  1.  ABSOLUTE  COORDINATE  SYSTEM 


KMC!  Of  DISIOHTiON  SCR  LENS  ON  TOTAL  PRESSURE  RiSE  OF  ROTOR  SET  NO. 


DISTORTION  SCREEN  DISTORTION  SCREEN 

NO  DISTORTION  CONFIGURATION  A.  <2  LOBEI  CONFIGURATION  B.  (4  LOBE) 


t 


M ' : ’ i ' I * 

; INLET  DYNAMIC  PRESSURE,  ^ , 


: f i "!  r : • • i • • • 

INCEPTION  * 


mm 


fyvyvyyWVY} 

• : : : : i • : ; : ; • • 

TIME  SCALE  CHANGE  — -► 


, TOTAL  PRESSURE  DROP 
i j-  ACROSS  SCREEN,  APy 

FlSBiffiTFi : ! i i 


a)  DOWNSTREAM  TOTAL  PRESSURE  PROBE,  PT  . IN  FREE  ST 


INLET  DYNAMIC  PRESSURE,  ^ 

-* 1 - 


N.  , ' 

rr — ' 


s-'s,  /*'v 


INCEPTION- 


1 MCOOd 


TIME  SCALE  CHANGE 


TttiHHWittHlW 


TOTAL  PRESSURE  DROP 
ACROSS  SCREEN.  APT  ' 


r r ir  :* 

il  1 -i 

I i ' i i 


. • _ 

* i r l 


\ : .1 


r j r i 

, i / i 


iK 
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Figure  25 

INCEPTION  OF  ROTATING  STALL  WITH  DECREASING 
INLET  DYNAMIC  PRESSURE 

ROTATING  DISTORTION,  ROTOR  RPM  = 1000 

DISTORTION  ROTATING  IN  SAME  DIRECTION  AS  ROTOR 
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INCEPTION  OF  ROTATING  STALL  WITH  DECREASING 
INLET  DYNAMIC  PRESSURE 
ROTATING  DISTORTION,  ROTOR  RPM  - 1000 
DISTORTION  ROTATING  IN  SAME  DIRECTION  AS 
ROTOR 
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Figure  26  (Cont.) 


TOTAL  PRESSURE  DROP  ACROSS  ROTATING  DISTOR- 
TION SCREEN  FOR  CONSTANT  MEAN  INLET 
DYNAMIC  PRESSURE,  qQ  -s  1 IN.  H20,  ROTOR 
RPM  - 1000 
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Figure  27 

TOTAL  PRESSURE  DROP  ACROSS  ROTATING  DISTORTION 
SCREEN  FOR  CONSTANT  MEAN  INLET  DYNAMIC 
PRESSURE,  q0  ^ 1 IN.  H20,  ROTOR  RPM  - 1000 
DISTORTION  ROTATING  IN  SAME  DIRECTION  AS  ROTOR 
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NOTE:  COORDINATE  SYSTEM 
FIXED  TO  BLADES 


Figure  29  FINITE  THICKNESS  ACTUATOR  MODEL  FOR  ANALYSIS  OF  DISTORTION 


OVERALL  TOTAL  PRESSURE  LOSS  COEFFICIENT. 


Figure  33  FLOW  TURNING  PERFORMANCE  EXTRAPOLATIONS,  STATOR  SET  NO  1 - 
VARIOUS  STATOR  STAGGER  ANGLES,  6SM 


OVERALL  TOTAL  PRESSURE  LOSS  COEFFICIENT,  AC 


INLET  SWIRL  ANGLE,  0vdeg 


Figure  34  OVERALL  TOTAL  PRESSURE  LOSS  EXTRAPOLATIONS,  STATOR  SET  NO.  1 - 
VARIOUS  STATOR  STAGGER  ANGLES,  8 cm 


ARROW  INDICATES  EXPERIMENTAL 
INCEPTION  OF  ROTATING  STALL 
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Fir  re  37  COMPRESSIBILITY  EFFECT  ON  STABILITY  CHARACTERISTICS  FOR 
STATOR  SET  NO.  4 WITH  h cm  = 28.2° 
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